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INTRODUCTION 


The  projects  described  here  fall  under  four  broad  subject 
areas:  thermal  leasing,  defect  properties,  damage  studies,  and 
surface  absorption  measurements.  Portions  of  this  work  have 
been  reported  at  the  Fifth  Annual  Conference  on  Infrared  Laser 
Window  Materials,  Las  Vegas,  Dec.  197S.  In  some  cases,  technical 
details  included  in  those  conference  proceedings  are  not  re- 
peated here.  Nearly  all  the  projects  reported  are  related  to 
work  performed  under  a previous  contract  (F19628- 7S-C- 0080)  , and 
the  technical  and  final  reports  of  that  program  include  some 
background  material  on  the  projects  described  below. 
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A Quantitative  Experimental  Investigation 
of  Thermal  Leasing 

C.  P.  Christensen,  W.  H.  Steier,  R.  Joiner 
Department  of  Electrical  Engineering 
University  of  Southern  California 
Los  Angeles,  California  90007 

Abstract 

A quantitative  comparison  between  experimental  measurements  and 
leasing  theory  in  CdTe  shows  good  agreement  in  both  first  order  leasing  and 
lens  abber ration. 


Thermal  leasing  and  optical  distortion  in  infrared  window  materials 
have  been  treated  theoretically  by  several  investigators^,  and  qualitative 
agreement  with  theoretical  predictions  has  been  demonstrated  experimentally 
in  some  cases.  However,  a detailed,  quantitative  comparison  of  theory  and 
experiment  has  not  previously  been  reported  and  is  necessary  if  composite 
window  structures  are  to  be  designed  to  compensate  for  thermal  leasing  aber  - 
rations. 

In  this  investigation  thermal  leasing  experiments  are  conducted  using 
s well  characterized  50  watt  CO,  laser  as  the  optical  source  and  relatively 
high  loss  (fi'^0.  2 cm  ) CdTe  as  the  leasing  medium.  The  laser  operates  in 
a fundamental  TEMqq  Gaussian  mode.  Although  the  available  laser  power  is 
not  sufficient  to  induce  measurable  lensing  in  low  loss  materials,  operation 
at  the  50  watt  level  is  desirable  since  parasitic  distortion  in  beam -handling 
optics  can  be  minimized,  laser  mode  control  is  facilitated,  and  optical  meas- 
urements are  simplified  at  reduced  power  levels.  The  emphasis  here  is  to 
obtain  an  experimental  check  on  the  theory  of  lensing  and  hence  well  character- 
ized experimental  conditions  are  stressed.  The  assumption  is  that  the  theory 
can  be  scaled  to  low  loss  materials  and  high  laser  powers  which  are  of  prac- 
tical interest. 

The  lensing  medium  is  uncompensated  polycrystalline  CdTe,  5 cm  dia,  by 
1.  3cnj  thick;  the  relatively  high  loss  is  caused  by  free  electrons.  We  assume 


that  these  free  carriers  do  not  effectively  change  the  material  constants 
such  as  K,  the  heat  conductivity;  Cv,  the  heat  capacity;  and  n,  the  index 
of  refraction.  We  have  calculated  the  effect  of  the  free  carriers  on  dn/dT 
and  found  the  effect  negligible. 

The  experimental  setup  is  shown  in  Fig.  1 and  uses  a rotating  mirror 
to  scan  the  beam  past  a pinhole  detector.  The  electronics  shown  enable  us 
to  select  any  scan  up  to  20  seconds  after  the  sample  is  first  illuminated. 
Typical  results  are  shown  in  Fig.  2 for  an  incident  laser  power,  P,  of  25  W; 
for  a sample  to  beam  waist  distance,  L,  of  1 m;  and  for  a detector  to  sample 
distance,  d,  of  1 m.  Fig.  2a  is  on  a slow  time  scale  and  shows  many  scans 
to  illustrate  the  build-up  and  decay  of  the  on-axis  intensity  as  the  leasing 
develops.  The  individual  scans  clearly  show  the  growth  of  the  side  lobes 
caused  by  the  aberrations  in  the  induced  lens. 

To  compare  these  results  with  theory  requires  knowledge  of  the  input 
Gaussian  beam  parameters.  Figure  3 shows  a slow  scan  of  the  input  beam. 
From  these  scans  taken  at  two  distances  from  the  laser,  the  position  of  the 
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beam  waist  can  be  determined.  Our  res\ilts  show  the  waist  to  be  96  cm 
from  the  flat  output  mirror  of  the  laser  and  a waist  spot  size  radius  of  2.  4 
mm.  The  theory  for  an  empty  optical  cavity  using  the  mirrors  shown  would 
predict  the  waist  to  be  at  the  output  mirror;  the  difference  emphasizes  the 
need  for  measurement  of  the  beam  parameters  and  is  likely  caused  by  leas- 
ing effects  in  the  discharge  tube.  The  beam  parameters  were  found  to  be 
independent  of  the  laser  power  output  over  the  range  of  powers  used  in  the 
experiments. 

The  calculations  of  the  temperature  profiles  caused  by  the  absorbed 

laser  power  were  done  by  evaluating  an  exponential  integral  which  is  appro - 

4 

priate  for  an  infinite  media  and  were  also  done  by  a difference  equation 
technique  which  allows  arbitrary  radial  boundary  conditions.  The  solutions 
agreed  well  at  times  short  compared  to  the  thermal  diffusion  time  to  the 
sample  boundary.  Only  radial  heat  flow  was  considered,  since  heat  flow  la- 
terally off  the  sample  surface  was  expected  to  be  less  than  1%  of  the  radial 
heat  flow. 
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The  iacrease  in  absorption  coefficient  with  temperatures  due  to  the  increase 
in  free  electron  density  was  included  in  the  temperature  calculations  and  its 
effect  on  the  infrared  transmission  of  the  sample  was  included  in  the  diffrac- 
tion calculation.  The  infrared  absorption  increases  at  about  3%  per  degree 
centigrade  above  room  temperature. 

The  beam  diffraction  calculations  used  a scalar  wave  theory  and  as- 
stimed  radial  symmetry  and  a Gaussian  input  beam.  The  linear  expansion 
of  the  sample  with  temperature  was  included  and  made  an  approximate  10% 
contribution  to  the  lensing.  The  stress-optic  effects  were  neglected,  since 
the  sample  was  polycrystalline,  and  this  effect  is  expected  to  be  small  in 
CdTe.  No  resonant  interference  of  the  beams  reflected  from  the  sample 
surfaces^  were  considered  since  in  the  experiment  the  sample  surfaces 
were  not  normal  to  the  beam  axis.  We  did,  however,  include  the  absorbed 
power  from  the  first  reflected  beam  inside  the  sample. 

Figure  4 shows  a typical  comparison  between  the  measurements  and 
the  theoretical  calculations  of  the  transmitted  beam  profiles.  The  theory 
is  surprisingly  good  in  predicting  beam  shapes  but  somewhat  poorer  in  pre- 
dicting absolute  intensities,  particularly  at  longer  times.  The  agreement 
was  in  general  this  close  with  all  of  the  data  taken.  The  theory  was  found 
to  be  fairly  sensitive  to  the  values  of  material  parameters  used.  In  this 
comparison  we  used  the  best  available  values  for  CdTe  (K  = 0.  07  W/cm  °C, 

= 1.  22  joule/cm^  °C,  dn/dT  = 1.  2 x 10"'^  / °C). 

In  Figures  5,  6,  7,  and  8 are  plotted  the  on-axis  intensity  time 
and  the  side  lobe  intensity  vs.  time  for  several  experimental  configurations. 
The  theory  fairly  accurately  predicts  the  on-axis  intensity  and  correctly 
predicts  the  general  trend  of  the  side  lobe  intensity. 

^ a)  M.  Sparks,  J.  Appl.  Phys.  4^  5029  (1971) 

b)  J.  R.  Jasperse  and  P.  D.  Gianino,  J.  Appl.  Phys.,  1686  (1972) 

c)  J.  Marburger  and  M.  Flannery,  Conf.  on  High  Power  Laser  Window 
Materials.  AFCRL  - 71  - 0592,  Special  Report  No.  127  (1971) 

^ a)  li.  H.  Skolnik,  B.  Bendow,  and  E.  F.  Cross,  Appl.  Optics  ^ 726 
(1974) 

b)  B.  Bendow,  L.  H.  Skolnik,  E.  F.  Cross,  Appl.  Optics  13,  729  (1974) 

^ H.  Kogelnik  and  T.  Li,  Appl.  Optics  1550  (1966) 

^ J.  Gordon.  K.  Leite,  R.  Moore,  S.  Porto,  and  J.  Whinnery,  J.  Appl. 
Phys.  3^  ^ (1965) 

^ B.  Bendow,  L.  H.  Skolnik,  E.  F.  Cross,  Appl.  Optics  ^ 729  (1974) 
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SLOW  SCAN  OF  LASER  MODE 


2 sec 


AXIS  INTENSITY 
P = 34  watts 


Measurement  of  Strain -Optic  Coefficients 
By  Acousto -Optic  Interactions  in  Window  Materials 
R.  Joiner,  W.  H.  Steier,  and  C.  P.  Christensen 
Department  of  Electrical  Engineering 
University  of  Southern  California 
Los  Angeles*  California 

Abstract 

The  three  independent  strain-optic  coefficients  (Pj^j*  Pj^2  P^^)  were 
measured  for  CdTe  and  GaAs  at  10.  6^  , using  an  acousto-optic  technique. 
The  results  for  GaAs,  excepting  Pj2»  S’gree  with  the  theoretical  predictions 
of  Bendow  et  al. 


In  order  to  evaluate  the  performance  of  a high  power  laser  window,  it 

is  necessary  to  know  the  strain-optic  coefficients  (p..)  of  the  window  mater- 

1 2 3 4 ^^ 

ials.  Recent  experimental  work  at  10.  6/i  ' ' has  resulted  in  measure- 

ments of  ^tress -optic  coefficients  (^-j)  and  stress -birefringence  coefficients 
- qj2))  • Both  are  related  to  the  strain-optic  coefficients 
b>  the  stiffness  tensor  (p^^^  = This  paper  presents  a third  method 

which  was  originally  devised  to  measure  strain-optic  coefficients  in  candi- 
date materials  for  acousto-optic  modulators.  Using  this  technique,  we  were 

able  to  determine  all  three  independent  coefficients  for  both  GaAs  and  CdTe. 

5 6 

The  experimental  technique  was  first  demonstrated  by  Dixon  ' at 

7 8 

laser  wavelengths  of  . 63^  and  1.  IS^ti  and  later  extended  to  10.  Gfx  . Fig- 
ure 1 shows  the  experimental  set-up  for  10.  6^  . A 2 Msec  pulse  of  high  fre- 
quency sound  travels  roundtrip  through  a fused  quartz  modulator,  a salol 
bond,  and  the  window  material,  deflecting  a fraction  of  the  HeNe  laser  beam 
and  the  CO^  laser  beam  on  each  pass.  These  deflections  are  maximized 

when  the  Bragg  condition  X 

sine^a  = 2A  (1) 

9 ^ 

is  satisfied.  6^  is  the  angle  formed  by  the  optical  and  acoustic  propagation 
vectors.  X and  A are  the  optical  and  acoustic  wavelengths  respectively. 

The  per  cent  deflected^®  is  given  by 


= sin . 


M,  p 


acj 


(2) 


L and  H are  the  length  and  height  of  the  acoustic  beam,  and  p^^  is  the  acou- 
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Stic  power.  The  acousto-optic  figure  of  merit,  is  a function  of  the  in- 

dex of  refraction,  n ; the  density,  p ; the  acoustic  velocity,  V ; and  the  ap- 
propriate strain-optic  coefficient,  p. 


n 

P 


(3) 


Note  that  in  equation  (2)  a longer  wavelength  implies  smaller  deflections  for 
a given  figure  of  merit. 

Simply  solving  equation  (2)  for  M-2  reqxiires  a knowledge  of  the  acous- 
tic power,  which  is  a function  of  acoustic  losses  in  the  material  and  the  qual- 
ity of  the  salol  bond.  By  assuming  small  deflections  and  taking  ratios  of  the 
deflected  pulses  shown  in  Figure  2,  a simple  formula  relating  to  a known 
figure  of  merit  for  quartz  can  be  derived  and  does  not  require  a knowledge 
of  the  acoustic  power. 

^2  / 10.6/i  (4) 

M2  (quartz)  \lj  I3  / \0.  633^  / 


For  GaAs  and  CdTe,  there  are  only  three  independent  strain-optic  co- 
efficients, pjj,  pj2*  P44.  For  longitudinal  sound  propagating  along  a 
<100>  direction,  we  can  measure  p^  by  using  light  polarized  parallel  to  the 
plane  containing  the  light  and  sound  propagation  vectors,  and  we  can  measxire 
by  using  a light  polarization  perpendicular  to  that  plane.  For  sound  pro- 
pagating along  a ^111>  direction,  the  two  polarizations  yield 


Pll  " ^ (Pll  + ^Pl2  ^P4*i^ 

I 

P12  = ^ <Pll  + 2 P12  - 2 P44> 


(5) 

(6) 


P44  calculated  either  by  subtracting  equation  (6)  from  equation  (5)  or 

by  substituting  p^  and  p^2  either  equation  (5)  or  (6).  Since  strain-optic 
coefficients  can  be  either  positive  or  negative,  a separate,  but  simple  ex- 

g 

periment  must  be  done  to  determine  relative  signs. 

Some  care  and  compromises  are  required  to  assure  that  the  experi- 
mental conditions  are  appropriate  with  the  assumptions  inherent  in  equation 
(4).  Equation  (4)  is  valid  for  small  deflections  and  for  acoustic  frequencies 


7 
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near  those  defined  by  the  Bragg  regime  . It  is  also  necessary  that  the  IR 
beam  divergence  be  less  than  the  acoustic  beam  divergence  . In  our  case, 
we  were  able  to  use  an  acoustic  frequency  and  to  adjust  the  IR  beam  size  to 
reasonably  meet  these  requirements.  However,  this  resulted  in  the  sound 
pulse's  interacting  at  any  one  time  with  only  a portion  of  the  IR  beam,  which 
produced  a deflected  IR  pulse  longer  and  of  lower  intensity  than  was  ex- 
pected. To  account  for  this,  a boxcar  integrator  was  used  to  measure  the 
integrated  deflected  intensity. 

Samples  of  CdTe  and  GaAs  oriented  <100>  and  <111>  were  measured, 
and  the  results  of  those  measurements  are  shown  in  Table  I.  The  GaAs  sann- 
ples  were  measvired  at  an  acoustic  frequency  of  150  MHz,  and  the  CdTe  was 
measured  at  50  MHz.  Figure  3 shows  a comparison  of  experimental  (cir- 
cles) and  theoretical  (solid  lines)  values  for  GaAs,  The  values  at  1.15m  were 
6 12 

measured  by  Dixon  . The  theory  of  Bendow  et  al.  extended  these  results 
to  longer  wavelengths  and  predicted  little  change  at  10.  6^  . Experimental 
results  concur,  with  the  exception  of 

In  an  attempt  to  measure  polycrystalline  ZnSe,  it  was  found  that  the 
grain  sizes  were  of  the  order  of  an  acoustic  wavelength  and  no  coherent  in- 
teraction between  light  and  sound  was  observed  past  a few  millimeters  into 
the  ZnSe. 

This  method  offers  some  advantages  over  birefringence  or  interfero- 
metric methods  in  that  i)  since  the  strain  is  an  acoustic  traveling  wave,  the 
difficulties  of  applying  a uniform  strain  to  the  sample  and  of  residual  strain 
in  the  sample  are  avoided;  ii)  the  stability  of  an  interferometer  is  not  re- 
quired; and,  iii)  all  three  coefficients  for  cubic  materials  are  measured. 
However,  the  signal  in  this  method  tends  to  small  since  it  is  proportional 
p while  in  the  interferometer  method  the  signal  is  proportional  to  p. 
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Theory  of  Thermally  Induced  Birefringence  in  Polycrystal- 
line Media 

M.  Flannery,  J.  H.  Marhurger 

During  the  initial  period  of  this  contract,  we  developed 
what  promises  to  be  an  accurate  theory  of  stress  induced  bire- 
fringence in  polycrystalline  media.  This  theory,  which  is  based 
on  work  of  Hershey  and  Dahlgren^  on  the  elastic  properties  of 

polycrystalline  media,  was  reported  at  the  Fifth  Annual  Con- 

2 

ference  on  IR  Window  Materials  and  also  published  as  an  Applied 
Physics  Letter.^ 

The  result  of  this  theory  is  an  expression  for  the  com- 
ponents of  the  effective  strain  optic  tensor  p^^  of  an  ideal 
polycrystalline  material  in  terms  of  the  single  crystal  values 
of  the  tensor  components  P—  and  of  the  elastic  coefficients 
c^j . The  theory  has  been  developed  only  for  materials  of  cubic 
symmetry  of  the  classes  43m,  432  and  m3m,  which  include  all  the 
currently  interesting  IR  window  materials.  For  these  classes, 
only  the  tensor  components  Pi2’  ^44  ^o^^vanishing . (We 

use  the  usual  6 dimensional  matrix  notation  for  the  symmetric 
tensors  main  results  of  the  theory  are: 

1.  The  combination  (p^^+2p^2)/^  is  equal  to  the  single 
crystal  expression  (p^ + 2p 2) elastic  coefficients 
also  obey  this  relation: 

(c*^  + 20^2^/^  " ^^11  2Cj2)/3  = bulk  modulus. 

2.  The  combination  p*^  *^12  satisfies 


where  p = - ^*2) 

3(3lt  + 4y)c^^ 

~ (9/t+  8u)p  + 6(K.+  2pyc^^ 

3.  The  effective  shear  modulus  p of  the  polycrystalline 
medium  is  the  only  positive  root  of  a quartic  equation  whose 
coefficients  depend  upon  c^2  ‘'44'  value  of  p 
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determined  from  this  equation  falls  between  the  Hashin- 

Shtrikman  bounds,  which  are  derived  from  very  general  considera- 
4 

tions . 

Having  determined  ^ , it  is  necessary  to  compare  the  per- 
formance of  a polycrys tal 1 inc  window  with  that  of  a single  crystal 
window.  This  is  not  easy  because,  even  for  cubic  materials,  an 
arbitrarily  oriented  window  does  not  possess  sufficient  sym- 
metry in  the  window  plane  to  allow  an  analytical  solution  of  the 
thermal  stress  problem.  Only  in  one  case  is  the  analysis  trac- 
table: When  the  window  is  parallel  to  a (111)  plane,  the  rela- 
tion between  stress  and  strain  in  the  plane  stress  approximation 
has  the  same  tensor  form  as  that  for  isotropic  media,  and  closed 
expressions  can  be  obtained  for  the  strain  in  a nonuniformly 
heated  window.  This  configuration  was  treated  by  F.  Horrigan 
(private  communication)  and  also  by  us,  and  the  results  were  re- 
ported in  the  Fourth  Annual  IR  Window  Conference  proceedings.^ 

The  results  for  a (111)  window  can  be  summarized  as  follows: 
The  induced  birefringence  effects  are  described  by  the  dif- 
ference between  the  phases  of  radial  and  azimuthal  polariza- 
tion components  of  the  incident  optical  field.  (We  are  assum- 
ing a circularly  symmetric  intensity  distribution,  centered  on 
a sel f- constrained  edge-cooled  window.)  In  the  short  pulse 
heating  regime,  where  thermal  conductivity  does  not  distort  the 
temperature  distribution,  the  phase  may  be  written  simply  as 

$ _ = F[g(r)  - 2f(r)] 

where  f(r)  and  g(r)  are  related  to  the  incident  beam  shape.  For 
a gaussian  beam  with  1/e  intensity  radius  a, 

f(r)  = 

g(r)  = [l  - exp(-r^/a^)]/ (r^/2a^) 

The  dimensionless  parameter  F is 


exp(-r^Za^) 


F = 


/3  lotWn-^K, 

2F 


2CJ1C44  + SkTc/j-c^P 
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P12  . 2p 
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- c 
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44 


•44 
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r 


where  8 = absorption  coefficient 

Iq  = axial  incident  intensity 
t = illumination  time 
L = window  thickness 
a = thermal  expansion  coefficient 
n = window  refractive  index 
1C  = bulk  modulus 
c'  = volume  specific  heat. 

In  the  isotropic  limit  c^^ 


'44 


12 


), 


P44 


P44  = 


Sj(p*^  - p*2)  , we  have 


pc 


= A 


( 9 ^Pfl 


12 


>1 


where  A is  the  factor  multiplying  the  curly  brackets.  Prior 
to  the  development  of  a polycrystalline  theory,  the  birefrin- 
gence contribution  was  estimated  by^ 


= A 


9 ^Pll  ■ Pl2^  I 

p (3  k + y)  j 

where  p^^^^  and  single  crystal  parameters. 

In  Table  I,  we  list  the  values  of  F/A,  Ppj./A,  and  F^/A. 

The  better  figure  should  no  longer  be  used  as  a figure  of  merit. 
The  others  are  proportional  to  figures  of  merit  for  (111) 
single  crystal  windows  and  for  polycrystalline  windows,  respec- 
tively. 

We  conclude  from  Table  I that  the  estimates  of  the  earlier 
crude  theory  are  closer  in  magnitude  to  the  polycrystalline  re- 
sults, although  the  predicted  sign  of  the  effect  is  frequently 
incorrect.  In  passing  from  (111)  single  crystal  to  polycrystal, 
the  magnitude  of  the  birefringence  is  reduced  by  up  to  a factor 
of  ten,  but  the  sign  is  not  changed.  In  the  important  cases  of 
Cap2  and  Bap2,  the  correct  polycrystalline  values  are  respec- 
tively factors  of  2 and  3 greater  than  the  previous  estimate, 
although  the  polycrystalline  forms  of  these  materials  still 
have  respectively  9 and  5 times  loss  birefringence  than  the 
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(Ill)  single  crystal. 


Table  I.  Thermally  induced  birefringence  parameters  for  win- 
dows fabricated  from  materials  with  cubic  symmetry  in  their 
crystalline  form. 


Material 

o 

> 

F/A 

pc 

NaCl 

- . 00532 

-0.0199 

-0.00346 

KCl 

. 00932 

-0.0668 

-0.00412 

KBr 

. 0104 

-0. 0684 

-0.00396 

CaF^ 

- . 0255 

0.420 

0.0470 

BaF, 

- . 0191 

0.286 

0.0555 

Ga? 

- . 00241 

-0. Cl  32 

-0.00387 

GaAs 

- . 00104 

-0. 0122 

-0.00365 

Ge 

. 00140 

0.0223 

0.00605 

Diamond 

- . 00375 

-0.0130 

-0.00272 

The  table  entries  are  proportional  to  the  phase  difference  in- 
troduced between  radial  and  azimuthal  polarization  components 
of  the  incident  optical  field.  Pq/A  is  the  approximation  used 
previously  to  estimate  the  magnitude  of  the  effect.  F/A  is  the 
exact  value  for  a (111)  single  crystal  window. 
theoretical  value  for  an  ideal  polycrystal  1 ine  window. 

Our  next  task  is  to  determine  whether  partial  orientation 
of  the  grains  can  reduce  the  birefringence  still  further. 

This  investigation  requires  a modification  of  the  entire  fiershey- 
Dahlgren  model  of  polycrystalline  materials,  and  is  now  in 
progress . 
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An  Electron  Microscope  Study  of  the  Structural 
Imperfections  in  II-VI  Compounds  (CdTe  and  ZnSe) 


1 


G.H.  Narayanan  and  S.H.  Rustomji 

The  goal  of  this  part  of  the  research  program  is  to  identify  and  charac- 
terize the  lattice  defects  and  other  structural  inhomogenieties  in  CdTe  and 
ZnSe  by  using  transmission  electron  microscopy  techniques.  The  work  in 
this  period  was  concentrated  on  the  study  of  P-doped  CdTe.  The  electrict.1 
and  optical  properties  of  this  material  have  been  also  investigated  in  our 
laboratories  by  Professor  Kroger's  group  and  by  Professor  Spitzer's  group 
respectively. 

It  is  known  that  the  native  defect  centers,  impurities  and  precipitates 
control  the  electrical  and  optical  properties  of  CdTe  and  that  the  precise 
nature  and  concentrations  of  these  defects  are  dependent  on  the  crystal 
growth  conditions  and  post -growth  annealing  treatments.  In  the  present  in- 
vestigation, therefore,  it  was  proposed  to  examine  the  defect  substructure 
of  CdTe:P  crystals  in  the  as-grown  condition  as  a function  of  the  P concen- 
trations and  further  to  investigate  how  the  defect  structure  was  influenced 
by  post-growth  annealing  treatments  under  well-defined  partial  pressures  of 
Cd  and  Te.  Also,  attempts  were  made  to  correlate  the  results  of  the  present 

studies  with  the  electrical  and  ootical  orooerties  measurements  carried  out  here 

at  use  on  materials  of  identical  chemical  composition  and  thermal  history. 

19  18  17 

CdTe  ingots  with  P concentrations  of  5 x 10  , 5 x 10  , 2 x 10  and 

16  -3 

5 X 10  cm  were  grown  from  melt  by  the  horizontal  Bridgman  technique  in 
our  laboratory  by  Dr.  B.V.  Dutt.  No  attempt  was  made  to  control  the  over- 
pressure during  crystal  gro'vth.  Samples  were  cut  from  each  ingot  and 
annealed  at  700°C  under  Cd  partial  pressures  of  2.4  x lo  1 x 10 


Hi 
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1 X 10  , 1 X 10  and  2 x 10  atms.  for  periods  of  time  long  enough  to 

achieve  equilibrium.  All  the  saiqDles  were  quenched  to  room  tenperature 
following  heat  treatments.  Both  as-grown  and  annealed  samples  were  first 
mechanically  ground  from  both  surfaces  to  produce  wafers  ~200v  in  thickness 
and  subsequently  cut  into  specimens  3 x 3 mm  in  size.  Final  thinning  to 
obtain  electron  transparent  specimens  was  done  by  a chemical  jet 
polishing  technique  described  previously  (1) . 

Fig.  1 (a)  through  (c)  shows  typical  bright  field  micrographs  of 

19 

CdTe:P  crystals  in  the  as-grcwn  condition  for  P concentrations  of  5 x 10  , 

18  16 

5 X 10  , eind  5 X 10  cm  respectively.  Evident  in  these  micrographs  is 

the  presence  of  fine  precipitates,  approximately  100-200  A in  size,  dis- 
tributed uniformly  throu^out  the  matrix.  The  density  of  the  precipitate 
is  higher  in  crystal  with  higher  P content.  Since  no  diffraction  effects 
attributable  to  the  precipitate  appeared  in  the  electron  diffraction 
patterns,  a positive  identification  of  the  nature  of  the  precipitate  phase 
could  not  be  made.  However,  the  apparent  dependence  of  the  density  of 
precipitates  on  P content  suggests  that  the  precipitates  might  be  some 
compound  involving  P.  There  are  two  known  phosphorous  compounds  in  the 
Cd-Te-P  system,  Cd^P^  and  worth  considering.  In  the  pre- 

sent case  the  precipitates  are  more  likely  to  be  those  of  Cd^P^  if  we  con- 
sider the  fact  that  P was  introduced  into  the  crystal  throui^  additions  of 

Cd^P^  into  the  melt  during  crystal  growth.  The  formation  of  these  compounds 

(2) 

have  been  speculated  also  by  Hall  and  Woodbury  in  their  investigations  of 
the  diffusion  and  solubility  of  P in  CdTe. 

/part  from  the  precipitates  the  as-grcwn  crystals  were  virtually  free 
of  other  types  of  defects  such  as  voids  which  result  from  the  condensation 


of  excess  vacancies  as  has  been  observed  in  some  vapous  grown  CdTe  . The 

4 2 

average  dislocation  density  in  all  the  crystals  studies  was  well  below  10/  /cm  . 


Defects  in  Annealed  Crystals 


Fig.  2 (a)  through  (d)  shows  typical  microstructures  of  CdTe  specimens 
19  -3 

with  [P]  = 5 X 10  cm  following  annealing  at  700°C  under  progressively 

-5  -3  -1 

higher  Cd  partial  pressures  -2.4  x 10  , 2 x 10  , 1 x 10  and  2 x 10 

respectively.  It  can  be  seen  that  at  the  lowest  Cd  pressure  (Te  partial 
pressure)  the  microstructure  reveals  the  presence  of  a moderately  high 
density  of  voids  distributed  uniformly  throughout  the  matrix.  A careful 
excunination  of  these  voids  shows  that  in  many  cases  they  are  partial Iv 
filled  with  a second  phase  (indicated  by  the  arrows) • These  voids  are 
polygonal  in  shape  and  range  in  size  from  700  to  1200  A.  Their  number  den- 
sity is  --S  X lO^Vcm^.  The  origin  of  these  partially  filled  voids  can  be 

( 4) 

understood  on  the  basis  of  the  mechanism  proposed  by  Shiozwa  et  al.  . It 
has  been  pointed  out  that  annealing  under  Te  over-pressures  (i.e.  very  low 
Cd  pressure)  could  lead  to  the  presence  of  equilibrium  concentration  of  Cd 
vacancies  that  are  in  excess  of  Te  vacancies  at  the  annealing  tenperature. 

On  cooling  to  room  tenperature , the  Cd  vaccincies  coalesce  into  microscopic 
voids,  with  the  Te  atoms  which  would  normally  occupy  the  lattice  sites  in 
that  region  of  the  crystal  partially  filling  this  void.  This  type  of  voids 
have  been  previously  observed  in  vapour  grcwn  CdTe^^^. 

-3 

In  specimens  annealed  under  a Cd  pressure  of  2 x 10  atms  there  vias  a 
slight  decrease  in  the  volume  faction  of  the  voids.  In  this  specimen, 
hcwever,  there  were  essentially  two  different  distributions  of  void  sizes  - 
a few  large  partially  filled  voids  of  size  500-1400  A coexisting  with  a 
fairly  uniform  distribution  of  large  numbers  of  extremely  fine  voids  100- 

O 

150  A in  size  were  found  to  exist  as  shown  in  Fig.  1 (b) . The  matrix  re- 
gions adjacent  to  the  larger  voids  were  found  to  be  virtually  free  of  the 
smaller  voids.  In  many  cases  the  larger  voids  were  found  to  be  preferentially 
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nucleated  on  matrix  dislocations.  The  origin  of  these  two  distinctly  dif- 
ferent size  distributions  of  voids  is  not  clearly  understood  at  the  pre- 
sent time.  In  addition  to  the  voids,  a high  density  of  precipitates  can 
also  be  seen  in  this  microstructure.  Annealing  under  still  higher  partial 
pressures  of  Cd,  did  not  result  in  the  formation  of  any  identifiable  voids. 

However  a distribution  of  fine  precipitates  was  observed  in  specimens 

-1  -2 

annealed  at  Cd  partial  pressure  of  1 x 10  atm.  but  not  at  1 x 10  . An 

example  of  this  precipitation  is  shown  in  Fig.  2.  (c)  . At  the  highest  Cd 

partial  pressure  used  in  this  study,  namely  0.2  atms,  the  formation  of  a 

few  voids  were  again  observed.  These  voids,  in  general,  were  found  to  be 

tetrahedral  in  shape  and  nucleated  predominantly  on  the  matrix  dislocations. 

They  are  believed  to  be  due  to  the  condensation  of  vacancies  of  tellurium. 

In  addition  to  Cd  partial  pressure,  the  total  concentration  of  [P]  in 

the  crystal  was  also  found  to  have  a pronounced  effect  on  the  s\ibstr\icture 

of  annealed  crystals.  When  annealing  was  carried  out  at  very  lew  Cd  partial 

pressure  viz;  2.4  x 10  ^ atm.,  the  volume  fraction  of  the  void-precipitate 

conplex  was  found  to  decrease  progressively  with  decreases  in  P content 

as  shewn  in  Fig.  3 (a)  and  (b) . At  a p , of  1 x 10  ^ atms  voids  were  still 

Cd 

observed  in  all  crystals,  except  the  one  with  the  lowest  P content,  namely 
16  3 

10  /cm  . However  the  void  volume  fraction  was  consistently  lower  than 

those  observed  in  crystals  under  low  p_..  At  intermediate  partial  pressure 

Ca 

i.e.  1 X 10  ^ atms,  precipitates  100-300  A in  size  were  observed  in 
crystals  containing  2 x 10^^  and  6 x 10^^  P/cm^.  Crystals  with  hi^er  P 
concentrations  were  virtually  free  of  any  defects  following  this  heat 


treatment.  Fig.  4 (a)  and  (b)  show  typical  microstructures  of  crystals 

with  2 X 10^^  and  6 x 10^^  cm  ^ of  P respectively  following  annealing 

-2 

under  a p^^  = 1 x 10  atms.  Annealing  at  higher  Cd  pressures,  0.1  atms; 
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produced  materials  vdiich  exhibited  little  or  no  precipitation  when  P con- 

18  17  3 

centration  was  either  2 x 10  or  2 x 10  /cm  . While  small  amounts  of  pre- 
cipitated phase  was  observed  in  crystals  with  6 x 10^^  P/cm^  and  5 x 10^^  P/cm^. 

The  results  of  the  present  electron  microscopy  studies  are,  in  general, 

(5) 

consxstent  with  the  defect  model  proposed  by  Selim  and  Kroger  based  on 

their  resiilts  of  Hall  effect  measurements  and  previous  diffusion  and  solu- 

(2) 

bility  measurements  of  Hall  and  Woodbury  . Using  the  data  presented  by 

Selim  and  Kroger  it  is  possible  to  construct  defect  isotherms  valid  foi 

700°C,  in  which  the  concentrations  of  the  defects  dominant  in  the  neutrality 

ccxiditions  are  plotted  as  a function  of  the  P concentration  for  various  Cd 

over  pressure.  The  effects  of  p and  total  P content  of  the  crystal  on  the 

Cq 

annealed  microstructures  can  be  satisfactorily  ejgjlained  on  the  basis  of 

these  isotherms.  At  very  low  Cd  partial  pressures,  i.e. : 2.4  x 10  ^ atms 

the  isotherm  predicts  Cd-vacancies  as  the  predominant  defect,  the  concen- 

18 

tration  of  which  is  estimated  to  be  of  the  order  of  10  /cm'  for  crystals 
19 

with  5 X 10  /cm  of  P.  When  these  crystals  are  cooled  to  room  teitperature , 

coalescence  of  the  excess  Cd  vacancies  can  occur  leading  to  formation  of 

voids,  half  filled  with  Te  atoms  which  would  normally  occupy  the  lattice  site 

in  that  region.  At  a p^^  of  10  ^ atms  the  isotherms  predict  a much  lower 

concentration  at  all  P concentration,  thus  accounting  for  the  decrease 

in  void  volume  fraction.  The  isotherm  also  shows  that  when  total  concen- 
18  3 

tration  exceeds  10  /cm  , P^  becomes  a dominant  defect  species,  their  con- 

18  3 19  3 

centrations  being  of  the  order  of  10  /cm  in  crystals  with  10  P/cm  . Precipitation 

of  these  interstitial  P can  occur  during  cooling  to  room  temperature,  which  would 

explain  the  presence  of  precipitates  in  crystals  with  10^®  and  10^^/  P/cm^.  For 

-2 

equal  to  or  greater  than  10  atms  V'  is  no  longer  a dominant  defect  species  at 


20 


all  P concentrations.  Consequently  no  voids  are  to  be  expected,  in  agree- 
ment with  experimental  observations.  However  the  isotherm  shows  at  this 

16 

p„,,  P.  to  be  the  dominant  defect  in  crystals  with  P contents  of  6 x 10 
Cd  1 

17  3 18  19  3 

and  2 X 10  /cm  . For  P concentrations  of  10  and  10  /cm  , P.  was  found 

1 

to  form  neutral  complexes  with  Based  on  these  isotherms  one  can  pre- 

dict that  the  precipitation  of  P^  might  occur  in  crystals  with  10^^  and 
17 

10  [P]  and  little  or  no  precipitation  in  crystals  with  higher  P content. 

This  is  in  agreement  with  the  electron  microscopy  observation.  When 

annealed  under  a p„,  of  10  ^ atms,  the  defect  isotherms  show  Cd?*  as  a 
Cd  1 

18  3 

dominant  defect  species  in  crystals  with  P content  up  to  ~10  /cm  ; beyond 

* 

this,  the  formation  of  complexes  are  predicted.  According  to 

this  model  therefore,  precipitates  are  predicted  in  the  annealed  crystals 
at  lower  P concentrations  which  is  what  was  observed  in  the  microstructure. 
It  should  be  pointed  out  that  distinction  between  precipitates  of  P^  cuid 
Cd**  are  diffucult  under  the  present  condition  owing  to  the  absence  of 
any  detectable  diffraction  effects  due  to  these  precipitate  phases  in  the 
electron  diffraction  patterns. 

During  the  remainder  of  the  contract  period  we  propose  to  carry  out 
the  following  studies ; 


(i)  an  investigation  of  the  microstructures  of  CdTe  doped  with  Ge 
both  in  the  as-grown  and  annealed  conditions. 

(ii)  an  investigation  of  the  microstructural  characteristics  of 
ZnSe,  doped  with  various  concentrations  of  Al  following  various 
annealing  treatments. 

The  electrical  and  optical  properties  measurements  on  these  materials 
have  been  carried  out  here  at  USC  by  Professor  Kroger  and  his  students. 

An  atteirpt  will  be,  therefore,  made  to  correlate  the  microstructural 
features  with  the  meaisured  physical  properties. 


1 


1 
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2x10  P/cm' 


700 C Annealed  CdTe:P;  p =2-4x10  atm 


700*C  Annealed  CdTetP;  p. 
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study  of  Defects  in  II-VI  Compounds 
F.A.  KrSqer,  F.A.  Selim  eind  A.  Ray 


CdTeiP 

The  work  on  CdTetP  is  completed.  A defect  model  is  proposed 
contcdning  P^^  and  P^  as  the  acceptors  at  high  and  medium  respec- 
tively, and  P^^  as  a triple  donor  at  low  Solubility  data  by 

Woodbury  and  Hall  are  explained  assuming  major  phos- 

phorus species  at  low  P^^f  the  acceptor  species  mentioned  above  as  the 
major  P species  at  high  P^^j*  Local  mode  studies  by  Dutt  and  Spitzer 
indicate  the  presence  of  both  P^^  and  ^ t>ut  not  of  P^.  Publication: 

F.A.  Selim  and  F.A.  Kroger  "The  Defect  Structure  of  P-Doped  CdTe",  sub- 
mitted to  J.  Electrochemical  Society. 


CdTe  + Ge 

Measurements  on  the  carrier  concentration  and  10.6ym  absorp- 
17  -3 

tion  of  CdTe  +10  Ge  cm  quenched  after  equilibration  at  700“C  in 

-2  -1 

various  cadmium  pressures  previously  showed  a..  ^ = 3 x 10  cm  at 

10.6 

-3  , , ,„-3  -1  ^ ,„15  -3  ^ ,,,13  -3 

c.  = 10  cm  and  a,  „ = 3 x 10  cm  at  10  cm  > c,  > 10  cm  . 

h 10.6  h 

Repetition  of  these  experiments  led  to  values  of  * 4 X 10  ^ cm  ^ 


after  annealing  at  medium  p„„. 

Cd 

We  have  grown  new  ingots  with  10 ^^Ge  cm  10 ^^Ge  cm  and 
(lO^^Ge  + 10^^In)cm  ^ and  are  ready  to  anneal  these  and  measure  10.6ym 


absorption. 


ZnSe 

Measurements  of  the  high-temperature  Hall  effect  and  resis- 
tivity of  ZnSe  doped  with  3 ppm  Al,  30  ppm  A1  and  300  ppm  Al,  equilibrated 
with  various  partial  pressures  of  zinc  show  oin  electron  concentration 
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F 


1/2  1/2 

c “ P™  [Al]  at  medium  p . The  samples  with  3 and  30  ppm  Al  show 
0 Zn 

a constant  electron  concentration  with  c » c , at  high  p 

e Al  ^ ^Zn 

After  cooling  to  room  temperature,  the  sanples  annealed  at 

800®C  at  high  and  medium  p show  the  same  electron  concentrations 

Zn 

-4 

found  at  high  temperature  with  a sharp  cut-off  at  p^^  = 3.15  x 10  atm 

for  the  Scunple  with  3 ppm  Al  and  = 8 x 10  ^atm  for  the  sample  with 

300  ppm  Al.  Below  the  cut-off  the  sanples  have  a high  resistivity.  The 

cut-off  is  probably  due  to  the  trapping  of  electrons  at  V'  . The  pro- 

Zn 

1/2  1/2 

portionality  of  c to  p [Al]  at  low  p is  explained  by  a defect 
0 Zn  Zn 

model  in  which  the  major  charged  defects  are  Al*  ,v"  or  Se" : 

Zn  Zn  1 


Zn  + V„  ->■  Zn„  + 2e' 
Zn  Zn 


K. 


ZnV 


[e'l  = K p [v"  ] •<  1/2K  D [Al*  ] 
Znv  ^Zn^  Zn  ' Zn  \TZn^  Zn 


The  tenperature  dependence  of  the  electron  concentration  at  constant  p, 
and  [Al] , 


Zn 


[e']  “ exp  + 0.5  eVAT 


shows  that  1/2H„  „ = -0.5  eV  or  H . = -1  eV 

Zn,V  Zn,V 


A sanple  doped  with  50  ppm  As  (an  acceptor)  annealed  at  800,  900,  1000°C 


under  various  p^^  shows  still  a negative  Hall  effect  at  high  and  medium 

p_  with  c “ but  now  the  temperature  dependence  is  in  the  opposite 

zn  0 Zn 

direction : 


28 


[e*] 


exp  -0.8  eV/kT 


OC 


indicating  a positive  energy  of  activation.  This  is  to  be  expected  for 
a defect  model  dominated  by 


[As-  ] = 2[V*-]  or  2[Zn.  ] 
Se  Se  1 


or  a similar  expression  with  instead  of 

Thus  H (or  H)  + H , = 1.6  eV  or,  with  = -1  eV,  H^(or  H„)  = 2.6  eV 

S F Zn,V  Zn,V  S F 

The  absolute  value  of  [e']  measured  in  the  sample  doped  with 
As  are  similar  to  those  reported  by  Smith  for  undoped  ZnSe.  On  the 
otlier  hand  the  activation  energy  found  by  Smith^^ , 1.6  eV  + 0.2  cV  is 
the  same  as  the  one  found  by  us.  Therefore  it  is  likely  that  Smith's 
crystal  was  unintentionally  doped  with  acceptors  at  a concentration 
not  too  different  from  the  As  concentration  in  our  crystal. 

ZnSe ;Ga 

Attempts  to  measure  the  high-temperature  Hall  effect  on 

samples  doped  with  3 cind  30  ppm  Ga  in  equilibrium  with  atmospheres  of 

known  p_  were  unsuccessful. 

Zn 

Crystals  doped  with  100  ppm  Ga,  quenched  after  annealing  at 

800°C  at  high  p„  were  n-type  with  a concentration  of  electrons  < con- 
Zn 

centration  of  Ga.  Crystals  annealed  at  low  p^^  had  a high  resistance. 

2) 

Robinson  and  Kun  found  p-type  conductivity  for  ZnSerGa  after  annealing 
under  p^^  indicating  acceptor  behavior  of  Ga.  On  the  other  hand  Tyagi 
and  Arora^^  report  n-type  behavior  (2  x 10^^  e'  cm  for  ZnSe  annealed 
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with  Ga+Zn  at  950-700°C.  Yet  acceptor  behavior  of  Ga  is  indicated 
also  in  this  case  by  the  fact  that  ZaSe+Ga  annealed  in  zinc  vapor  has 

4) 

a lower  electron  concentration  than  pure  ZnSe  under  tlie  same  conditions 
It  is  likely  that  Ga  shows  amphoteric  behavior,  acting  completely  or 
partly  as  an  acceptor  at  high  p^,^  but  ar;  a donor  at  medium  and  low 
Further  work  is  needed  to  clarify  the  defect  model. 
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Ahsorpt  ion  Studies  of  CdTo  and  ZiiSe 
W.  G.  Spitzer  and  B.  V.  Putt 

This  section  deals  with  absorption  mechanisms  in  CdTe  and 
ZnSe,  particularly  impurity- related  absorption  and  free  carrier 
absorption. 

A.  Background  and  Summary  of  Our  Previous  Work. 

Studies  of  infrared  absorption  of  doped  CdTe  and  ZnSe  have 
been  the  subject  of  our  previous  progress  reports.  The  dopants 
chosen  were  A1  and  P,  as  they  are  light  in  mass  compared  to  the 
host  atoms.  Tliis  choice  enables  identification  of  the  defects 
responsible  for  the  observed  absorption,  since  the  free  carrier 
absorption  is  reduced  due  to  the  prevalent  sel f - compensat ion 
observed  in  II-VI  compounds.  Free  carrier  absorptions  can  be 
introduced  by  annealing  in  controlled  partial  pressures  of  the 
constituents.  Throughout  the  investigations  the  primary  long 
range  goal  has  been  to  identify  the  various  absorption  mech- 
anisms and  their  spectral  range  of  influence.  It  is  generally 
known  that  the  various  factors  limiting  the  transparency  of 
the  II-VI  compounds  are  multiphonon  absorption,  absorption  due 
to  precipitates  or  clusters  of  impurities  or  native  defects,  and 
inter-  and  intra-band  carrier  absorption  involving  electronic 
transitions.  .All  these  factors,  except  the  mu  1 1 i ph<^non  pro- 
cesses, are  largely  influenced  by  the  conditions  under  wliich 
the  crystals  are  grown  and/or  post-growth  ' annea  Is . A substan- 
tial effort  has  been  directed  to  preparing  our  own  crystals  and 
studying  the  influence  of  post-growth  anneals  on  the  contribu- 
tions of  the  various  absorption  mechanisms.  The  following  sec- 
tions briefly  summarize  our  work  on  CdTe;Al,  ZnSeiAl,  and  CdTe:P 
previously  reported  in  detail. 

fl)  CdT e : A1 . It  is  usual  to  add  donor  elements  In  or  Cl  to 
CdTe  to  promote  se 1 f- compensat ion . This  results  in  a high  re- 
sistivity material  that  is  quite  transparent  in  the  infrared. 
Various  defect  models  involving  the  presence  of  unassociated 
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^”cd’  ^Cd’  ‘Associates  ^ “^h'e'''Cd^ 

have  been  proposed.  These  models  were  largely  based  on  indirect 
evidence  obtained  from  transport  measurements.  The  infrared 
localired  vibrational  mode  technique  is  a more  direct  way  of 
identifying  the  defect  spec  ies,  a Ithough  it  is  limited  to  ele- 
ments considerably  lighter  than  the  host  crystal  atoms.  There- 
fore, the  choice  for  a donor  is  essentially  restricted  to  the 
element  Al.  The  electrical  behavior  of  A1  should  be  similar  to 
In  in  CdTe.  We  therefore  grew  CdTc.-Al  crystals  which  were  in- 
tentionally doped  with  an  [All  = lO^^cm  . The  introduction  of 
counterdopants  such  as  Cu,  Ag  and  Au,  usually  by  diffusion,  was 
found  to  be  useful  for  the  identification  of  the  different  Al- 


related  defect  species'.  Full  details  of  this  work,  giving  pos- 
sible defect  assignments,  has  been  reported  elsewhere.  When 
CdTe.'Al  was  annealed  in  Cd-rich  conditions,  free  carrier  ab- 
sorption was  observed.  The  absorption  is  nearly  temperature 
independent,  and  shows  a(v)  'v  v and  v ^ after  appropriate  an- 
neals. The  V behavior  was  explained  in  terms  of  polar  op- 
tical mode  scattering  and  the  v ^ behavior  was  accounted  for  in 
terms  of  ionized  impurity  scattering.  There  were  two  major 
concerns  that  were  noted:  the  absence  of  temperature  dependence 

in  the  a 'v  v ^ measurements  and  the  better  than  expected  results 

18  - ^ 

that  the  theory  gave  at  the  higher  carrier  densities  10  "cm 
where  plasmon -phonon  coupling  is  expected  to  cause  the  theory 
to  be  invalid. 


(2)  CdTe : P . This  work  was  started  primarily  to  understand  the 
self-compensation  mechanisms  through  the  use  of  LVM  tecimiques 
in  acceptor-doped  II-VI  compounds.  The  materia]  was  grown  in 
our  laboratory  by  using  the  horizontal  Bridgman  method.  Pre- 
vious progress  reports  contain  details  of  our  absorption  mea- 
surements on  this  system.  The  work  is  now  complete  and  has 


been  submitted  for  publication. 


The  major  findings  arc: 


(a)  there  are  two  P defect  centers,  one  P„,  or  P.  center  having 

1 c 1 

T_j  symmetry  and  the  other  a P-comp)c.x  with  axial  symmetry; 
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(M  the  P,^.^  or  is  a shallow  acceptor,  while  the  P-complex 
appears  to  be  a neutral  center;  (ci  the  shallow  acceptors  P,j,^ 
or  P^  are  present  in  unassociated  form  and  native  donors,  such 
as  Cd  interstitials,  account  for  the  observed  self-compensation. 

(3)  ZnSe : A1 . The  crystals  used  in  this  work  were  from  Czochralski- 
grown  ingots  from  Eagle-Picher  Co.  A set  of  ZnSetAl  was  an- 
nealed in  Zn-rich  conditions  and  some  free  carrier  absorption 
I measurements  were  previously  reported.  The  free  carrier  ab- 

sorption shows  a(v)  ny  V behavior  and  is  nearly  temperature 
independent.  These  results  are  now  analyzed  and  will  be  pre- 
sented in  a later  section. 

We  also  carried  out  LVM  measurements  on  ZnSe:Al  with  Cu, 

Ag  or  Au.  Infrared  LVM  bands  attributed  to  pair  defects  of 
Al-Cu,  Al-Ag  or  Al-Au  were  observed.  These  were  very  similar 
to  the  corresponding  defects  in  CdTe.  The  results  have  been 
published  elscw'here . 

B . Work  Done  During  the  Present  Progress  Report  Period. 

(I)  CdTe : P . The  work  on  this  system  has  been  completed  and 
presented  for  publication. Some  of  the  findings  were 
already  mentioned  in  the  above  background  and  summary  section. 

In  order  to  better  understand  some  of  the  unrff^iolved  details  of 
the  sources  of  defect  - induced  absorption  in  CdTe:P,  we  have 
initiated  studies  of  the  Zn  Cd, _ Te:P  and  Mg  Cd,  Te;P  binary 
. alloy  systems.  The  crystals  were  all  grown  using  the  horizontal 

Bridgman  method.  The  absorption  spectra  are  obtained  at  So'^K 
on  the  as-grown  samples.  The  free  carrier  absorption  is  low 
and  did  not  cause  any  problems  of  measuring  the  LVM  bands  from 
P-related  defects.  Preliminary  indications  are  that  the  two 
systems  behave  differently  with  respect  to  the  originally  ob- 
served P defects.  The  results  are  being  analyzed  and  will  be 
I reported  later. 

i 

I 


[D  InSc : A 1 . In  continuation  of  our  previous  work,  we  have 
now  carried  out  detailed  calculations  of  free  carrier  ahsorp- 
tion  of  ZnSe  by  using  tlie  tlieorv  for  juilar  optical  mode  scatter 
ing.‘^  The  spectral  dependence  of  the  atisorption,  the  nearly 
linear  dependence  on  carrier  density,  tiie  independence  of  the 
A1  concentration  and  the  magnitude  of  tiie  absorption  cross 
section  are  all  in  accord  with  this  theory.  However,  as  noted 
earlier,  the  free  carrier  absorption  measurements  showed  nearly 
temperature  independent  absorption  between  3D0°K'  and  8()^K.  This 
fact  cannot  be  reconciled  with  the  present  theoretical  model 
for  polar  optical  mode  scattering  absorption.  This  behavior  of 
annealed  ZnSe:.Al  is  very  similar  to  that  noticed  earlier  for 
CdTe:Al.^  Full  details  of  our  work  are  in  preparation  for  pub- 

f* 

licaticn. 
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High  Power  Infrared  Laser  Damage  Studies  , 

M.  Bass 

This  part  of  the  program  is  concerned  with  CO^,  DF  and  HF 
pulsed  laser  induced  damage  studies.  We  have  shown  that  the  surface  j 

and  bulk  damage  threshold  of  several  materials  at  10.6  jam  depends  i 

not  only  on  such  beam  parameters  as  pulse  duration  and  irradiation  j 

area,  but  also  on  tlie  sequence  of  pulses  used.  These  results  are 
described  in  Ref.  1 (see  Attachment  A)  and  point  up  the  necessity  for 
re-examination  of  the  conclusions  drawn  about  damage  mechanisms 
from  measurements  of  relative  damage  thresholds. 

Having  obtained  reliable  performance  from  the  infrared  modula- 
tor, (see  Fig.  1),  we  are  starting  to  secure  the  data  necessary  for 
survival  curve  analyses  of  the  10.6  )jm  damage  mechanisms  in  IR 

window  materials.  This  powerful  diagnostic  technique  was  described 

2 3 

originally  by  Milam  £t  al.  ’ In  order  to  obtain  more  accurate  data 
recording,  a computer  controlled  data  acquisition  system  was  requested 
for  the  current  program.  As  a result  of  delays  in  contract  negotiations, 
this  system  was  not  ordered  until  April  of  1976  and  is  not  expected  to  be 
available  for  use  before  mid-July.  It  appears,  therefore,  that  the  sur- 
vival curve  studies  will  be  only  partially  completed  by  the  end  of  the 
current  project. 

In  this  program,  we  have  started  the  first  pulsed  laser  damage 
studies  in  the  2-5  |im  region  of  the  spectrum.  This  part  of  the  project 
has  been  conducted  jointly  with  the  Naval  Weapons  Center  (NWC)  at 
China  Lake,  California.  Our  initial  efforts  were  aimed  at  operating  and 
characterizing  the  HF/DF  TEA  laser  at  NWC.  We  found  that  the  laser's 
energy  output  was  so  irreproducible  on  a pulse  to  pulse  basis  that  con- 
stant monitoring  was  required.  The  dependence  of  the  pulse  waveform 
on  the  wavelength  of  oscillation  made  it  difficult  to  obtain  reliable  data 
and  showed  the  need  for  single  line  operation.  In  addition,  a believe-  ' 

able  measurement  of  the  focal  spot  diameter  has  not  yet  been  obtained  ^ 
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hooauso  tho  linht  burns  tiu'  pinholos  used  to  scan  tin*  beam.  llowt'\er, 
under  tbesi'  eircumsi.inees  vvi'  lia\e  be<'n  abli-  to  study  the  ri'lali\c 
damage  thresholds  of  sexeral  materials  of  interest.  Our  results  for 
this  initial  effort  are  listed  in  Table  I with  several  other  possibly 
related  material  properties.  Table  II  lists  the  scattering  data  obtained 
for  2-5  (Jnii  window  materials  in  this  work.  From  the  fact  that  several 
samples  have  a poorly  defined  threshold  in  a 1 on  1 experiment,  we  can 
say  that  inclusion  damage  appears  to  be  prevalent.  It  is  also  clear  that 
current  finishing  techniques  do  not  produce  damage  resistant  surfaces. 

In  fact,  the  differences  between  bulk  and  surface  thresholds  are  in 
most  cases  greater  than  would  be  predicted  on  the  basis  of  field  enhance- 

4 

ment  at  polishing  imperfections.  If  "vee"  grooves  were  present  all 
over  the  surface,  field  enhancement  could  explain  the  damage  data  but 
microscopy  does  not  show  a sufficient  density  of  such  grooves.  The 
fact  that  field  enhancement  can  not  adequately  explain  the  difference  in 
damage  thresholds  may  indicate  that  the  surfaces  are  contaminated  by 
absorptive  polishing  compounds  or  that  as  a result  of  polishing  they  are 
highly  absorptive.  Damage  studies  of  etched  surfaces  and  calorimetric 
measurements  of  surface  and  bulk  absorption  are  planned. 

To  obtain  more  quantitative  data,  the  laser  is  being  more  accu- 
rately characterized  and,  if  possible,  single  line  operation  with  dam- 
aging intensities  will  be  obtained.  Several  other  candidate  materials 
will  be  studied  and  correlations  sought  between  2-5  Hm  damage  thresh- 

I 

olds  and  other  material  properties.  • j 

We  have  attempted  to  secure  the  most  up  to  date  and  promising  ' 

materials  for  damage  testing.  We  have  obtained  the  mose  recent  CVD 
grown  7nSe  from  Raytheon,  RAP  processed  KCl  and  NaCl  from  Hughes 
Research  Laboratories  and  pressed  forged  KCl  from  Honeywell.  The 
latter  was  found  to  be  recrystallized  but  will  be  studied  while  a new  i 

sample  is  prepared.  Initial  data  on  the  new  7nSe  was  presented  in  j 

Ref.  5.  We  are  currently  preparing  this  material  using  an  etching  j 

procedure  developed  at  Hughes  and  survival  curve  studies  are  in  j 


ah 


progress.  The  RAP  processed  samples  are  being  finished  and  a RAP 
processed  KBr  sample  is  being  sent  to  us  from  the  Naval  Research 
l.aboratories.  We  have  arranged  with  the  Hughes  Research  Hab  1o 
provide  us  with  samples  of  improved  1 0.  h jam  AH  coatings  for  damage 
resistance  testing.  These  are  expected  to  start  arriving  in  May. 
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Table  I.  Relative  Damage  Thresholds  Using  "All- line"  DF  "TEA"  Laser  and  Other  Properties 


r.edian  values  cf  widely  separated  minimum  fluxes  for  oarage  and  maximum  for  no  damag 
that ''these  media  are  inclusion  damage  limited. 


Table  II.  Scattering*  in  2-5ym  Window  Materials 
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waveform  before 
modulator 


Waveform  trans- 
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ATTACHEMENT  A 


I!  II-  JOliKNAl  Ot-  Mil  IKOMtS,  1>BKUAKV  1970 


The  Dependence  of  the  Pulsed  10.6-Min  Laser  Damage 
Threshold  on  the  Marnier  in  Which  a Sample  is  Irradiated 


M BASS  AN11  K.  M LLUNG 


AO^rrael-  fhp  [hieshold  tor  pulsed  10.6-uiti  lasei-induced  surface  and 
bulS  damage  in  iianspateiil  nuleiials  has  been  shown  to  depend  on  the 
manner  in  which  the  samples  arc  irradiated.  When  a .site  is  inadiated 
several  times,  starting  at  a very  kiw  intensity  and  increasing  the  flux  .a 
few  percent  on  successive  shots  until  damage  occurs,  the  threshold  is 
often  higher  than  if  each  site  were  inadiated  only  once.  In  addition, 
it  has  been  shown  that  the  dauiagc  threshold  also  depends  on  the  sam- 
ple tested  and  so  the  use  of  relative  damage-threshold  measurements 
to  identify  the  damage  mechanum  should  be  reexamined. 


Thresholds  tor  pulsed  .aser  damage  to  transparent  media 
have  been  determined  in  u \ariety  of  expenments  The  impor- 
tance oi  the  beam's  spatial  mode  and  the  pulse’s  temporal 
wave''  rm  have  been  discussed  (1)  Tl'.cre  has  also  been  ex- 
terisi.e  discussion  of  the  relation  bti.vc  'ii  the  size  of  the  ir- 
radiated volume  and  ihv  measured  th'cshoid  12'  In  this 
correspondence,  we  poini  out  the  impoilance  ol  the  num- 
ber of  times  a site  is  irradiated  in  determining  the  rcpo.ted 
damage  threshold  at  10.6  urn  and  now  this  fact  aftects  com- 
parisons between  material'. 

We  have  repotted  experm’-’it , in  which  each  site  was  irradi- 
ated only  on^e  and  the  intensity  was  m-  reased  at  each  succes- 
sive sue  until  damage  was  produced  with  some  degree  of  con- 
sistency (l  e .'0  percent  i ' ihe  sites  are  damaged!  [.^|.  Allen 
f/  a!,  reported  damage  t.'^ooholds  obtained  w cn  inch  .site  wa.s 
irradiated  .several  tunes  iiarting  at  a very  low  intensity  and  in- 
creasing the  flux  a tew  pcr  'int  on  successive  st;ols  until  dam- 
age occurred  14| . Since  our  one  shot  on  one  site  ( 1-on-l)  and 
then  several  shots  on  one  site  t,V-on-l)  expenments  gave  dif- 
ferent results  for  the  same  materials,  we  have  undertaken  a 
direct  comparison  of  tl..'se  procedures  on  the  same  samples. 

Our  experiment  is  described  m detail  in  ld|  and  is  summa- 
rized m Table  I The  results  of  out  1-on  1 and  /V-on-1  com- 
parison are  given  in  Table  II.  Clearly,  when  the  diameter  of 
the  irradiated  region  is  ~I0(I  iim  the  10.6-pm  laser  damage 
threshold  of  /nSe  is  independent  ol  the  manner  of  ii radiation 
Kt  I and  NaCI,  on  the  other  hand,  have  higher  damage  thresh- 
olds when  ircasureil  in  an  V on  I expeiimriit  than  in  a l-on  i 
experiment  This  increased  ihresholil  is  observed  lor  '-.‘>0  pei- 
cent  ol  the  sites  studied  in  both  surface  and  bulk  damage  ex- 
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periments.  No  dependence  on  interpulse  inter  . . -s  found 

for  a range  from  5 s to  >1  min. 

As  in  1 .1) , we  also  monitored  the  transmitted  pulse  waveform 
distortion  during  these  experiments  We  found  that  in  an  A'-on-l 
expeiiment  no  distortion  could  be  detected  Uurmg  any  pulse 
prior  to  the  one  which  produced  ,i  spark  and  left  re.sidual  damage 
in  the  sample.  The  distortion  of  the  last  pulse  was,  in  all  re- 
spects, similar  to  that  seen  in  a 1-on-l  damage  measurement. 

The  incrcas,'  in  surface  damage  threshold  in  an  /V-on-l  ex- 
periment might  be  due  to  eithei  a cleaning  of  the  surface  by 
the  laser  at  low  intensities  (5)  or  a smoothing  out  of  the  sharp 
edges  of  polishing  scratches  by  some  form  of  “lasei  polishing.” 
If  sharp  edges  are  removed  and  “vee”-like  scratches  become 
more  nearly  cylindrical,  Bloembergen  (61  would  predict  a ra- 
tio of  damage  intensities  for  /V-on-I  to  I -on- 1 of 

- 

4 

where  (■  - is  the  optical  frequency  dieleciric  constant,  and 
R is  the  ratio  ot  the  square  of  Bloembergen's  field-enhance- 
ment factors  for  “vee"  grooves  to  that  for  cylindrical  grooves 
Thus,  for  NaCt  where  ii  = 1 49,  we  should  have 
and  for  KCl  where  n-  14b  we  should  have  /C ice,  = 2.45. 
From  Table  II  we  see  that  the  experimental  ratios  arc  .t  06' 
3.23  and  2.36.  respectively  in  other  words,  the  change  in 
Ihe  surface  damage  threshold  between  the  A'-on  1 and  1-on-l 
experiments  does  not  disagree  with  the  "laser-polishing”  ex- 
planation The  tact  that  Ihe  surface  damage  threshold  of  ZnSe 
IS  independent  of  the  type  of  experiment  may  be  due  to  the 
fact  that  Ihe  damage  mechanism  in  this  material  is  absorption 
by  particulate  inclusions  13)  This  mechanism  is  not  sensitive 
to  local  field  enhancement. 

The  foregoing  model  might  explain  the  results  for  surface 
damage,  but  does  not  suffice  when  wc  consider  the  results  for 
bulk  damage.  The  increased  damage  flux  in  an  A'-on-l  experi- 
ment for  hulk  damage  is  both  unexpected  and  troublesome. 
Consider,  for  example,  the  fact  that  the  ranking  of  the  two 
NaCI  samples  reverses  when  one  goes  from  a 1-on-l  threshold 
to  an  A!-on-l  value  Consider,  in  addition,  comparing  kCl  lv> 
NaCI.  For  each  type  of  experiment  one  gets  a different  rela 
live  damage  threshold.  The  fact  that  Ihe  relative  damage 
threshold  also  depends  on  which  samples  arc  compared  sug- 
gests that  the  use  ol  such  comparisons  to  identify  the  damage 
mechanisms  (7|  he  reexamined. 

In  studying  the  statistics  of  laser  damage,  at  I 06  and  0 69 
pm  the  site  wa.s  irradiated  several  times  with  pulses  of  fixed 
intensity  |8|  The  irradiation  continued  until  eilliei  damage 
occurred  or  until  some  predeleimined  number  ol  pulses  had 
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been  tried.  \i  any  particular  intensity  level  the  distribution 
of  ine  nunibcr  of  pulses  required  to  produce  damage  was  given 
by  the  binomial  distribution  [S),  [9],  This  is  a characteristic 
of  a process  having  a con«'ai'.  probability  of  occurrence  on 
each  trial.  If  low-level  irrada^tion  caused  an  increase  in  the 
damage  threshold  ol  the  man  als  studied  at  1 .Oo  and  0.69jum. 
then  the  damage  probabiliij  would  not  have  been  constant 
tTom  pulse  to  pulse.  The  nrol'abihty  for  damage  would  have 
decreased  on  each  siiec 'Rsn  ,•  pulse  so  tlial  if  damage  dul 
not  occur  on  ihe  first  puls-,  a would  veiv  likely  not  have  oc- 
curred at  all,  Ihe  process  would  have  appearcii  ihresholdlike 
tnstead  of  statistical.  W.-  feel,  therefore,  that  the  increased 
threshold  m an  A'-on  1 experiment  at  10.6  ^n\  does  not  con- 
tlici  with  the  statistical  properties  of  1,06-  and  0.69-ktm  laser 
oamage.  In  fact,  the  statistics  observed  at  the  shorter  wave- 
lengths suggest  that  there  would  be  no  increa.se  in  threshold  in 
an  ,V-on-l  experiment  performed  with  either  an  Nd  or  a ruby 
laser. 

The  measured  10.6-/;m  laser  damage  threshold  for  both  bulk 
and  surface  damage  in  some  materials  has  been  shown  to  de- 
pend on  the  procedure  used  I his  study  has  also  shown  that 
the  relative  ranking  of  materials  according  to  their  damage 
thresholds  depends  on  whicn  -.amples  arc  compared  aiui  on  the 
experimental  procedure.  \s  j result,  the  use  of  such  rankings 
to  identify  the  laser  damage  mechanism  should  be  reexamined. 

Note  Added  m Proof:  No  increase  in  bulk  damage  threshold 
for  the  sample  of  NaCi  (llaishaw  Spectrophotometer  Grade) 
was  observed  in  an  A'-on- 1 experiment  as  compared  to  a 1-on-l 
experiment  performed  with  a ruby  laser 

.AC  KNOW  !,H)GMK.NT 

The  authors  wish  to  thank  C.  Tang  and  M Chang  for  their 
assistance  m this  work.  They  also  wish  to  thank  V.  Wang,  who 
participated  in  the  initial  experiments. 

Rl r PREMCKS 

[1|  ,M  Bass,  "Laser  induced  damage  to  nonlinear  crystalline  materials," 
NB.S  Special  Publ.  341,  pp.  90-96,  Ifec.  I97ii. 

|2|  L.  G.  De.Shagcr,  B.  E.  Newnam,  and  K M Leung.  "The  role  of 
coating  defects  in  laser-induced  damage  to  thin  films,"  NBS  Spe- 
cial PubL  387,  pp.  114-123,  Dec  1973 
|3|  K M.  Leung,  .M  Bass,  and  A.  G.  J halhm  Villavcide.  '"Damage  to 
l<).6  pm  window  materials  due  to  <"Oj  TEA  laser  pulses,"  pie- 
senled  al  the  1975  laiser  Dam.ige  Svmp  . Boulder,  f t),  lulv  I97s, 
and  to  appear  in  the  proeeedings  ol  that  meeling., 

|4)  S.  D.  Allen,  M.  Braiinstcm.  (. . Giuliano,  ami  V Wang  '"Pulsed  ("Oj 
laser  damage  studies  of  RAP  grown  KCI,"  NB.S  Special  Publ  414, 


pp.  66-75,  Dec.  1974. 

(5)  T.  L.  Barber.  “Degradation  of  laser  optical  surfaces,”  ELO.M  Kep. 
5237,  Mar.  1969. 

(6)  N.  Bloembergen,  “The  role  of  cracks,  pores  and  absorbmg  imlu 
sions  on  laser  induced  damage  threshold  at  surfaces  of  tianspateni 
dielectrics,"  Appi.  Opr.,  vol.  1 2.  pp,  661 -664,  Apr.  1973. 

[7j  E.  Y.iblonovitch,  "Optical  dielectric  strength  of  alkali-halide  crys 
tals  obtained  by  laser  induced  breakdow  n," /I pp/  P/ii's.  /.err  ,vol. 
19,  pp.  495-497,  Dec.  1971. 

|8i  M.  Bass  and  H.  H.  Barrett.  "Avalanche  breakdown  and  the  probahi 
listic  rijiuie  oflaser-induced  damage," /EE/V  Quar’liini  hleilron  . 
voL01-8.pp,  338-.34  3,Mar.  1972. 

[9|  M.  Bass  ami  1).  W.  E'tadin,  "Surface  and  bulk  lasetslamage  statistics 
and  the  identification  of  intrinsic  breakdown  processes,"  /E'E'E"  J. 
Quantum  Electron.,  vo\.  QE-9,  pp.  890-896,  Sept.  1973. 


Surl'iici'  Ah  so  rpr  i on  Snuiios  with  Atoustio  I'roho  Ti'ohn  i qiies 

.1.  II.  Parks  and  P.  A.  Rockwoll 

This  final  report  discusses  the  completion  and  accomplish- 
ments of  an  experimental  program  which  applied  acoustic  surface 
wave  technology  to  the  problem  of  measuring  the  weak  radiative 
absorption  occurring  within  a thin  layer  at  the  surface  of  a 
dielectric  solid.  Details  of  these  studies  including  analysis 
and  application  to  KCl  samples  are  presented  in  the  following 
jiages  in  the  form  of  two  technical  papers  presently  acceptc'd  for 
publication  in  the  .lournal  of  Applied  Physics. 

As  a general  overview  of  the  program,  the  surface  acoustic 
wave  techniques  were  found  to  be  an  extremelv  sensitive  tool 
which  could  be  used  to  unravel  the  surface  and  bulk  absorption 
of  10. 6y  radiation  in  prepared  KCl  materials.  It  was  observed 
that  surface  absorption  as  small  as  a 1 0 ^ could  be  directly 
measured  in  the  presence  of  rather  largo  hulk  absorption, 

Sg^^lO  ^cm  A theory  was  derived  to  extract  these  radiative 

parameters  from  experimental  data  and  in  fact  could  fit  the  data 
over  a wide  range  of  experimental  conditions.  The  experimental 
sens  i t ji  V i ty  , coupled  with  these  analytical  methods,  suggest  the 
capability  of  extending  measurements  of  the  surface  absorption 
to  0 £ 10  '.  Experimental  techniques  were  developed  which  allowed 
the  surface  waves  to  be  fluid  coupled  to  an  arbitrary  material 
under  study.  These  techniques  are  straightforward  and  should 
provide  the  flexibility  to  apply  the  measurement  technique  to  a 
wider  class  of  materials.  It  should  be  noted,  however,  that 
polycrystalline  materials  for  which  the  average  grain  size  is 
comparable  to  the  acoustic  wavelength  will  introduce  significant 
acoustic  wave  scattering.  In  this  regard  lower  acoustic  fre'- 
c^ucncies,  for  wliich  the  wavelength  is  comparaltle  to  SO-lOOnm, 
have  been  observed  to  propagate  on  polycrystalline  KCl.  The  ex- 
tension of  similar  measurements  to  ZnSe  appears  a realistic  pos- 
sibility. 

To  summarize,  the  principal  accomplishments  of  this  program 


I lu  I luU' : 


1.  Application  of  surface  acoustic  wave  Icchnology  to  the 
measurement  of  weak  surface  absorption. 

2.  Derivation  of  a tlieory  which  was  shown  to  extract  radi- 
ative absorption  parameters  directly  from  experimental 
data . 

3.  Development  of  techniques  to  couple  surface  waves  to 
arbitrary  materials  which  in  itself  broadens  the  ma- 
terials applicability  of  the  mcasureraent  technique. 

This  technology  is  presently  being  transferred  to  the  China 
Lake  Naval  Weapons  Center,  where,  under  the  guidance  of  Dr.  H. 
Bennett  and  Dr.  D.  White,  it  can  possibly  be  put  to  expanded 
use . 
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ABSTRACT 


The  application  of  an  Lcoubtic  technique  to  measure  surface 
radiative  absorption  in  highly  transparent  materials  is  presented. 
Dispersive  phase  changtiS  in  acoustic  surface  \^lave  piropagation  are 
induced  by  the  absorption  of  10.6  um  pulsed  laser  radiation  on  KCl 
surfaces.  Analysis  based  upon  an  absorbing  surface  layer  model 

provides  a jt  15%  determination  of  both  the  total  surface  absorption 

-4  -3  _i 

a - 4.  5 X 10  and  a bulk  absorption  coefficient  8 =2x10  cm  for 

13 

a mechanically  polished/chemically  etched  surface.  After  mechanical 

“4 

repolishing,  this  sample  yields  a = 11.3  X 10  . The  sensitivity  of  this 

technique  indicates  the  possibility  of  measuring  a surface  absorption  as 
-5 

small  as  cr  <10  and  also  obtaining  details  of  the  radiative  absorption 
depth  distribution. 

PACS  numbers:  78.  50.  - w,  81.80.  rx 
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The  weak  radiative  absorption  occurring  in  laser  window  materials 
represents  a practical  problem  limiting  the  application  of  high  power 
lasers.  Crystal  growth  procedures  h:ive  developed  the  capability  of 
growing  highly  pure  materials  with  bulk  absorption  coefficients  approaching 
intrinsic  levels.  As  a result,  major  research  efforts  have  been  directed 
Coward  studying  the  residual  surface  absorption  occurring  within  1-10  gim 
of  the  crystal  surface.  We  present  results^  of  a new  technique  to  measure 
absorption  of  CO^  laser  radiation  on  both  mechanically  and  chemically 
polished  KCl  samples.  These  results  indicate  (i)  high  sensitivity  to 
surface  effects,  (ii)  the  ability  to  separate  surface  <ind  bulk  absorption 
with  measurements  on  a single  sample,  .and  (iii)  the  possibility  of 
dete rrrvining  the  detailed  depth  dependence  of  the  absorptive  properties. 

The  theory  of  acoustic  surface  wave  detection  of  radiative  absorp- 

2,3 

tion  has  been  well  established  in  quantitative  studies  on  quartz  and 
3 

lithium  niobate.  The  fundamental  idea  is  that  the  temperature  change 

induced  in  the  material  by  the  radiative  absorption  results  in  a change  in 

the  velocity  of  a continuous  acoustic  surface  wave  propagating  through  the 

absorbing  region.  For  a fixed  acoustic  path  length  21,  along  the  x-axis, 

3 

the  velocity  change  Av  manifests  itself  as  a phase  change  given  by 

s 


L,  Av  (x,t) 

4tt  V ^ " s 1 

V Ii  J _ V J 

so  O so  A 


where  v is  the  frequency  and  v^^  is  the  unperturbed  velocity  of  the  surface 
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temperature  change. 

In  highly  transparent  materials  the  absorption  coefficient  0 within 

1-10  lam  of  the  surface  can  be  significantly  greater  than  that  in  the  bulK. 

The  extent  of  this  surface  absorption  depends  on  surface  preparation 

proceuures,  immediate  environmental  impurity  adsorption,  and  perhaps 

the  effect  of  surface  dislocations  induced  during  sample  preparation. 

The  simplest  model  that  embodies  these  effects  is  an  absorption  layer 

model  in  which  an  absorbing  surface  region  with  3 = 3^^  extends  a 

distance  d into  the  material,  at  which  de[)th  a transition  to  the  bulk  value 

3 - 3 occurs.  For  example,  in  KCl  at  the  10.  6 |am  laser  wavelength, 

Ij 

one  m igh  t expec t 3 , ~ f 0 ctTi  , d ~ 1 , and  3^^  < 1 0 cm  . This 

set  of  surface  parameters  results  in  a tot.il  surface  absiirption 
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cr  = = 10  , which  is  typical  for  mechanically  polished  surfaces. 

This  absorption  layer  model  is  used  to  calculate  the  temperature  change 
AT(r,z,t)  induced  by  radiative  absorption.  We  consider  a laser  pulse  of 
temporal  width  Xp,  gaussian  profile  of  width  (HWW)  a which  propagates 
along  the  z-axis  and  is  absorbed  by  the  material.  For  pulsewidths  con- 
sidered in  this  work,  the  temperature  variation  in  z,  characterized  by  the 
depth  dependence  of  B(z),  is  much  greater  than  radial  variation  diarac- 
terized  by  the  gaussian  width  a.  Radial  diffusion  is  then  negli^jible. 
Consequently,  the  radial  te.mperature  profile  is  identical  to  the  gaussian  laser 
pulse,  and  the  temperature  AT  is  found  by  soivir."  a cnc-dimensional  heat 
equation.  (Rigorous  solutions  to  the  three  dimensional  heat  equation 

justify  this  approximation.  ) The  material  is  considered  an  infinite  slab 

t 

of  thickness  b,  and  all  surface  heat  losses  are  neglected.  The  solution  is 

most  conveniently  written  in  series  form  AT(z,t)  = T^  (t;  0^  , d , 3 ^ ) c os 

where  k nn  /b. 
n 

The  dispersive  velocity  ciiange  [ Av  / v ] follows  from  a 

S S u K 

s 

4 

perturbative  treatment  of  the  acoustic  wave  e<|uation  in  which  the  elastic 
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coefficients  exhibit  spatial  variation  ihrouf^h  tlieir  weak  (—  10  ^/^'k) 

S 

tempo  la  til  ri-  ilepeiuleuce  . We  aj,suiue  lhal  lli.-rtnal  ilillusioii  ha.i  eMt-iule.l 

the  subsurface  temperature  li  is  I ri  l)ul  i i>n  to  a ileplh  '■  A for  which  a 
volume  perturbation  theory  is  appropriate.  Although  the  details 

of  this  calculation  are  not  essential  to  the  present  discussion,  it  is 
instructive  to  describe  the  functional  form  of  the  solution; 


^ Av  (x,  t) 
s 

L V J 
so  X 


E T (t;0  d,3  ) A(k  X ) e 

^ n s B ns 
n = 0 


X ^ 


J 


In  2 


(2) 


The  function  A approaches  a constant  for  small  values  of  its  argument; 

its  physical  significance  is  that  it  introduces  dispersion  for  terms  such 

that  ^ Experimental  measurements  of  this  dispersion  are  useful 

in  determining  the  absorption  [>arameters  since  each  X provides  a 

different  weighted  sum  of  the  T . The  coefficients  T determine  whether 

n n 

Av^/v^^  is  sensitive  to  and  d individually,  or  simply  their  product 
a - S^d,  as  follows.  Detailed  information  on  the  spatial  variation  of 
B(z)  is  contained  only  in  the  coeffici  c-nts  with  large  n.  For  tiires 
<■  10  usee,  Eq.(.^)  has  si  uni  f i i ,i  at  < ont  ri  but  i ons  fi-cim  siuh  terms,  ,md 
phase  change  irieasu  rement.s  will  determine  d and  ,B  indiqjendently.  I'or 
times  > 100  usee,  however,  most  of  the  high-n  terms  have  decayed, 
and  consequently  the  acoustic  phase  change  is  only  sensitive  to  cr. 

Experimental  phase  change  measurements  on  KCl  have  l>een 


.=10 


obtained  using  the  apparatus  shown  in  Fig.  1.  The  signal  lron\  a 


I 

i 


s 


I 


1 

I 


i 
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stabilized  VHF  generator,  amplified  to  1 watt,  launches  surface  waves 

cn  a quartz  transducer  crystal  using  iuterdigital  transducers.  The  surface 

7 

wave  is  fluid-coupled  onto  the  KCl  sample,  interacts  with  the  radiatively 
induced  temperature  distribution,  and  is  then  fluid-coupled  back  to  the 

transducer  crystal,  xhg  reliability  of  this  fluid-coupling  process  was 
previously  shown  in  Ref.  7 and  is  reaffirmed  by  the  present  results. 

The  transmitted  signal  enters  the  right  channel  of  the  double  balanced 
mixer.  The  left  channel  receives  a reference  signal  from  the  generator. 

The  variable  delay  line  is  adjusted  such  that  a transient  phase  change 

of  the  transmitted  signal  by  an  amount  Ao(t)  yields  a mixer  output 

signal  sin  A9(t),  where  is  typically  IbO  mV.  The  peax 

value  of  the  chase  change  is  measured  e.t  the  time  T , the  end  of  the  laser 

p’ 

pulse.  These  peak  values  are  detected  and  held  by  a peak  detector,  and 
averaged  by  a boxcar  integrator.  The  pulse  repetition  rate  is  chosen  to 
maintain  the  r.sme  average  laser  power  for  measurements  at  each  j-.ulsewidt''. 

T . 

p 

A (O2  laser  fiulsc  of  a,f)5  watts  is  focusoil  tiiroug.h  a liole  in  tl'.e  trans- 
ducer crystal  and  is  incident  on  a Harshaw  KCl  sample.  Acoustic  frequencies  of 

124.  b,  62.1,  and  2H.0  MHz  were  used  for  which  X -14.  1,28.2  and 

« ] 

62.  5 i_im  respectively.  Using  a focused  laser  beam  of  width  u ~ 0.  2 iiitm, 

3 

the  m easured  vivlues  of  f are  . 2 1 , . 2 b,  and  . 1 8 5 at  these  th  ree  ^ 

frequencies,  respectively.  Figure  2 shows  representative  oscilloscope 
trnctus  of  A9(h)  for  two  pulsewidths  at  121.  b MHz.  The  100  msec  [>eak 


r 


phase  change  in  Fig.  2a  is  60  mrad  corresponding  to  Av  /v  = -2.2  X 

s S o 
S 

per  watt  for  - 1.753  x 10  cm/sec.  The  peak  j)hase  change  of  2.1  mrad 

in  Fig.  2b  was  obtained  for  a 0.5  msec  55uiJ  laser  pulse,  and  is  about  teti 

times  greater  than  the  measured  rms  noise  of  '1^0.2  mrad.  This  establishes 

the  sensitivity  of  the  acoustic  technique  and  indicates  the  ability  to  study 

lower  absorption  coefficients  and,  if  necessary,  shorter  pulsewidths  than  are 

presented  here.  Experimental  uncertainties  arise  form  absolute  uncertainties 

in  the  measurement  of  the  following  quantities:  'v/±S%  each  in  the  laser  power, 

the  value  of  V^,  and  possibly  the  fluid  coupling  scheme,  and  'i-±io%  in  the  ex- 

3 

perimental  determination  of  fg.  These  lead  to  an  overall  absolute  uncertainty 
in  these  measurements  of  '^'±15%. 

The  data  set  taken  on  a sample  surface  which  was  both  mechanically 

Q 

polished  and  chemically  etched  is  shown  in  Fig.  3a  (referred  to  as  chemical 
polish).  The  curves  representing  the  theory  are  discussed  below. 

Following  these  measurements , the  same  KCl  sample  was  repolished  | 

mechanically  in  our  laboratory  using  standard  techniques,  and  a second  j 

set  of  measurements  were  taken.  These  mechanical  polish  results,  j 

i 

shown  in  Fig.  3b,  represent  the  two  lowe  r acoustic  frequencies^  tliis  case  ' 

fi,  is  .33  ami  .25  for  (i2.1  and  28.0MHz,  respectively.  Comparison  of  the  data  in  | 

U 

Figs.  3a  and  3b  indicates  that  the  chemical  etch  reduces  the  surface  absorption  o ; 


,! 
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by  a factor  of  ~3.  Similar  results  have  been  observed"  in  calorimetric 
measurements. 

A least  squares  computer  program  varied  the  two  parameters  of  the 

absorption  model,  a - 3^d  and  3^,  to  obtain  the  optimum  fit  of  Eqs.  (1) 

and  (2)  to  the  data.  As  is  clear  from  Fig.  3,  the  theory  represents  an 

adequate  description  of  the  data  except  for  the  shorter  times  and  lower 

frequencies.  It  is  believed  this  short  time  discrepancy  arises  from  a 

breakdown  of  the  volume  perturbation  calculation  of  Av  /v  when  the 

s so 

temperature  depth  dependence  becomes  comparable  to  the  acoustic 

penetration  depth  which  is  ~ 2.  2 for  KCl.  A separate  fit,  made  to  the 

124.  5 MHz  data  alone,  is  represented  by  the  curve  plotted  in  Fig.  3a. 

Because  of  the  shorter  penetration  depth  the  theory  is  adequate  at  that 

frequency  even  for  the  shortest  times.  The  data  for  both  lower  frequencies 

in  Fig.  3a  were  fit  simultaneously  using  pulsewidths  t > 10  msec.  The 

P - 

average  value  of  cr  based  on  these  two  independent  fits  to  the  chemical 
polish  data  is  a=  (4.  S'- . 5)  1 0''^  with  3 „ = 2 X 1 0'^cm  ' \ A simultaneous 

fit  was  also  done  lor  the  data  shown  in  Fig.  3b  using  the  same  [)uls ewidths , 
yielding  o ' 11.3  X 10  and  3j^  =2x10  cni  for  the  tnechariical  polist-i 
data.  These  results  for  a and  3^  are  consistent  with  calorimetric 
measurem  ents^^  on  these  and  similar  samples. 

Figure  4 shows  llie  reduced  phase  change  Acp/(vf|^)  versus 

laser  (julsewidth  on  the  i Iicm  ically  polisheil  surface.  Tl)e  reduced 
phase  change?  is  convenie?nt  in  discussing  dispersion  .since  A'p^^  idiminate?s 
the:  trivial  linear  fre:f)uency  depende»nce  in  Fq.  (1)  and  also  the  uninte*re>sting 


.S.3 


effects  arising  from  different  acoustic  profiles  at  each  frequency.  The 


divergence  of  the  data  for  shorter  times,  when  the  subsurface  temperature 
gradient  9AT/9x  becomes  comparable  to  AT  / X , indicates  the  presence  iif 
a surface  absorption  layer  even  after  significant  ir\aterial  (~10  gun)  w.is 
removed  during  the  etching  process.  Such  a temperature  gradient  could 
only  arise  from  a rapid  depth  variation  of  3(z)  on  the  scale  of  several 
microns.  The  importance  of  this  figure  is  to  establish  the  observation  of 
dispersion  well  outside  the  experimental  uncertainty  even  though  tlie 
present  theory  is  insufficient  to  describe  data  for  longer  wavelengths 
and  shorter  times. 

The  present  results  indicate  the  utility  and  sensitivity  of  this  acoustic 

technique.  Since  the  absorption  of  70  watt,  100  irisec  laser  [uilse  on 

-4 

a surface  with  cr  ~ 3 X 10  induced  a phase  change  of  Ap  - 60  m rad  at 

124.  5 MHz,  a readily  measureable  phase  change  of  1 tn  rad  would  result 

-3 

from  a surface  with  cr  ~ 10  . Dispersive  measurements  requiring  data 

_4 

at  lower  frequencies  are  less  sensitive  and  have  a lower  limit  of  o'  ~ 1 0 
based  on  the  28  MHz  results.  This  limit  is  comparable  to  the  lowest 
surface  absorption  measured  on  chemically  polished  samples. Both 
limits  could  be  further  reduced  by  o[ie  rating  ,al  higher  laser  power.  In 
order  Lo  measure  3 and  d independently,  it  is  aen*ssary  to  use  pulse- 
widlhs  T <10  psec.  Present  THA  lasers  (CO  , HF/Dl-')  provide  short 

p ~ 2 

[uilses  with  bufficic-nt  energy  for  such  measurements.  I'he  breakdown  vif 

the  [iresent  theory  for  short  times  rioes  not  .'(■jireseul  a practical  limitation 

4 

of  this  technique  since  theories  presently  exist  whicli  describe  surface 

S I 


wave  propagation  in  layered  media.  Such  a theory  could  be  used,  for 
example,  in  a study  of  radiative  absorption  in  films  of  thickness  < 1 pm 
if  laser  pulsewidths  < 1 psec  are  used  to  reduce  the  effect  of  thermal 


diffusion.  The  sensitivity  ol  this  technique  suggests  immediate  application 
to  several  important  measurements:  (i)  Details  of  B(z)  near  surfaces, 

(ii)  Impurity  diffusion  into  alkali-halide  surfaces,  and  (iii)  Absorption  in 
thin  dielectric  or  conducting  films. 

The  authors  wish  to  thank  S.  Allen  and  A.  Timber  of  Hughes 
Research  Laboratories  for  the  KCl  samples  and  their  preparation.  They 
also  thank  D.  White  and  R.  King  of  China  Lake  NWC  for  transducer 
fabrication,  and  D.  Smith  for  his  invaluable  assistance  in  writing  the 
computer  programs. 
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FIGURE  CAPTIONS 
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Fig.  2 


Fig.  3 


Fig.  4 


Ex|)i>  I'iuuMila  1 .ippa  iM  lud  tor  .sludyim;  rad  i.i  Li  vi-  al>si>rption 
at  I'iCl  surtaci;s  using  acoustic  techniques. 

Oscilloscope  traces  of  the  mixer  output  showing  the  tratisient 
phase  change  induced  on  chemically  polished  KCl  by  10.6  gm, 

70  watt  laser  pulses,  (a)  4 mv/div,  20  msec/div,  = 100  msec; 
(6)0,4  rav/div,  2 msec/div,  ip  = 0.5  msec.  Tlie  acoustic  f:e- 
quency  is  124.5  MHz. 

Phase  change  data  versus  laser  pulsewidth  at  various 
frequencies  for  (a)  chemical  polish  and  (b)  mechanical 
polish  KCI  surfaces.  Bars  indicate  J_  1 5%  experimental 
uncertainty.  Solid  curves  are  theoretical  fits  to  the  data. 

Phase  change  dispersion  versus  laser  pulsewidti  's 
indicated  by  the  increasing  divergence  of  as  t 

decreases.  These  data  were  taken  on  a chemically  polished 
surface;  experimental  uncertainty  is  _f  1 h'/o.  Dashed 
curves  are  shown  as  guides  to  the  i-yi'. 
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SURFACE  ACOUSTIC  WAVE  VELOCITY  DISPERSION 


INDUCED  BY  RADIATIVE  ABSORPTION* 
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ABSTRACT 


An  analysis  of  the  dispersive  propagation  of  surface  acoustic 
waves  through  a time  dependent  thermal  gradient  induced  by  radiative 
absorption  is  presented.  The  dispersive  phase  change  of  the  acoustic 
wave  is  shown  to  depend  on  the  details  of  the  surface  absorptive 
properties  of  the  material.  In  particular,  measurement  of  such  induced 
phase  changes  can  determine  the  parameters  which  characterize  the 
surface  depth  dependence  of  the  radiative  absorption  coefficient.  This 
analysis  is  applied  to  two  cases:  (i)  a strongly  absorbing  material  and 

(ii)  a highly  transparent  bulk  material  with  an  absorbing  surface  layer. 
The  relative  sensitivity  to  bulk  and  total  surface  absorption  is  discussed, 
as  well  as  the  radiative  pulse  conditions  for  which  the  absorption 
coefficient  depth  profile  can  be  determined.  This  technique  is  applied  to 
published  results  for  10.6  um  radiative  absorption  on  KCl  surfaces  for 
which  the  total  surface  layer  absorption  and  the  bulk  absorption  coefficient 
are  individually  measured.  Additional  applications  to  related  surface 
measurements  are  discussed. 


PACS  Numbers ; 78.  50.  - w,  81.80.  +x,  68.  90.  + y 
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I.  INTRODUCTION 

The  application  oi  surface  acoustic  wave  techniques  to  detect  and 

measure  surf.ice  r.idialive  a l>s  orj)t  i un  has  been  des»  i-il>etl  in  sev«Tal 

previous  pul>l  ii  a t i v>us  . The-  emphasis  in  (his  e.irliei"  work  ei-ulereci  on 

developing  sensitive  techniques  to  measure  radiatively  induced  acoustic 

phase  changes  on  arbitrary  materials,  and  deriving  a detailed 

4 

theoretical  description  of  the  thermal-surface  wave  interaction  . / .though 

the  observation  and  utility  of  induced  dispersive  propagation  was  discussed 
in  Refs.  3 and  4,  the  primary  concern  was  to  develop  a quantitative 
description  of  the  simplest  interaction,  which  occurs  when  the  tempera- 
ture is  uniform  over  the  depth  profile  of  the  surface  wave.  In  this  paper 
we  present  an  analysis  of  dispersive  surface  wave  propagation  arising  from 
the  time  dependent  thermal  gradient  induced  by  radiative  absorption. 

The  use  of  surface  acoustic  wave  velocity  dispersion  to  probe 
material  surface  properties  has  received  considerable  attention,  and  a 
review  of  these  applications  is  available  in  Ref.  5.  In  general,  the 
physical  basis  of  dispersion  arises  from  a perturbation  of  the  elastic 
coefficients  which  exhibits  an  appreciable  gradient  within  the  surface 
wave  penetration  depth.  The  energy  of  a surface  acoustic  wave  is 
largely  confined  to  within  a wavelength  of  the  surface,  so  that  probing 
the  material  surface  with  several  different  acoustic  frequencies  provides 
information  about  the  depth  variation  of  the  perturbation.  In  most  cases 
which  have  been  studied  the  perturbation  was  static,  and  an  accurate 
determination  of  the  gradient  relied  on  the  use  of  m.any  different 


b.’ 


frequencies.  This  can  be  impractical  when  it  is  desired  to  study  the 
surface  of  an  arbitrary  material  on  which  surface  wave  excitation 
introduces  additional  experimental  limitations  such  as  restricted 
frequency  range  or  surface  contamination.  In  contrast,  the  thermal 
gradient  resulting  from  radiative  absorption  evolves  in  time  through 
heat  diffusion.  In  this  case,  measurement  of  the  time  dependence  of 
dispersive  effects  provides  an  experimental  degree  of  freedom  which 
eliminates  the  need  for  many  frequencies  and  greatly  simplifies  data 
taking.  Since  the  process  of  heat  diffusion  near  the  surface  can  be 
described  reliably,  the  measurement  and  analysis  of  dispersive  surface 
wave  propagation  can  be  applied  to  unravel  the  depth  variation  of  the 
radiative  surface  properties  of  materials. 

This  paper  will  develop  the  analysis  necessary  for  the  application 
of  these  techniques.  In  Sec.  II,  the  acoustic  phase  change  is  derived  by 
relating  the  surface  wave  velocity  dispersion  to  the  temperature  distribu- 
tion arising  from  radiative  absorption.  Section  III  treats  two  limiting 
examples:  (A)  a strongly  absorbing  material  and  (B)  a highly  transparent 

material  with  an  absorbing  surface  layer.  The  theory  for  example  (B) 
is  applied  in  detail  to  the  problem  of  measuring  weak  surface  absorption 
on  KCl  samples.  Finally,  Sec.  IV  discusses  applications  to  several 
related  problems. 
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n.  DISPERSIVE  PHASE  CHANGE 

It  was  shown  in  Ref.  4 that  perturbed  surface  wave  propagation 
in  a region  heated  by  radiative  absorption  could  be  accurately  described 
by  the  thermally  induced  time  dependence  of  the  acoustic  phase  given  by 


^(t) 


2 Try  f 

V J 

so  path 


^v^(x,  t) 


v 

so 


dx 
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where  v is  the  frequency,  v^^  is  the  unperturbed  surface  wave  velocity, 
and  Av^(x.t)  is  the  change  in  the  surface  wave  velocity  due  to  the  thermal 
variation  of  the  elastic  coefficients.  The  integral  is  taken  over  the  total 
acoustic  path  length.  The  general  experimental  configuration  is  shown 
in  Fig.  1.  Laser  radiation  is  incident  on  the  surface  m.idway  along  the 
acoustic  path  of  length  2L.  The  laser  beam  direction  is  normal  to  the 
surface,  along  the  z-axis.  The  temperature  dependence  of  the  elastic 
coefficients  is  small  enough  10'^/°K)  that  it  is  appropriate  to 
linearize  the  dependence  of  the  acoustic  wave  vector  K = Zny/v  , ana 

subsequently  the  surface  wave  velocity,  on  the  induced  temperature 
change  AT  ( r,  t). 

In  general,  there  is  a contribution  to  (j>(t)  which  results  from  the 
change  in  total  acoustic  path  length  due  to  thermal  expansion.  In  the 
case  under  consideration,  absorbed  laser  radiation  will  heat  a local 
region,  but  the  surrounding  cooler  part  of  the  crystal  will  constrain  the 
thermal  expansion,  thereby  inducing  internal  stresses.  The  constrained 
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thermal  expansion  is  treated  by  replacing  the  mathematical  increment 

dx  by  a physical  distance  in  the  crystal  given  by  dx[l  +(9u/3x)].  The 

thermoelastic  strain  (Su/bx)  is  a measure  of  the  change  in  the  distance 

u between  surface  elements  due  to  the  nonuniform  thermal  stresses,  and 

can  be  calculated  under  the  assumption  of  plane  stress.  It  was  pointed 
4 

out  previously  that  the  constrained  thermal  expansion  present  within  a 
heated  region  having  dimensions  much  less  than  the  propagation  path 
generally  results  in  a contribution  to  the  change  in  ^(t)  of  less  than  a few 
percent.  This  path  length  expansion  will  be  omitted  in  this  paper.  In 
addition,  for  the  thermal  perturbations  under  consideration,  we  assume 
the  resulting  thermoelastic  stress  has  negligible  effect  on  the  wave  velocity. 

Under  these  conditions,  the  transient  phase  change  induced  by 
radiative  absorption  becomes 


so  o 


L r (x.t) 


assuming  an  axially  symmetric  temperature  distribution.  The  subscript 
that,  in  general,  the  integrand  depends  on  the  acoustic  wave- 
length. The  line  integral  is  taken  through  the  center  of  the  heated  region. 
The  factor  f^  accounts  for  the  incomplete  overlap  of  the  acoustic  beamwidth 
and  the  temperature  profile.  It  should  be  pointed  out  that  since  the  depth 
or  z dependence  is  of  central  importance,  we  consider  the  case  in  which 


9dT/dx  « AT/X  , because  surface  wave  reflections  from  the  perturbed 


II 

region  are  negligible  under  the  experimental  conditions.  In  the  derivation 
of  Eq.  (2),  it  if  assumed  that  the  time  dependence  of  the  perturbation  is 

4 

slow  compared  to  an  acoustic  wave  period.  Equation  (2)  was  derived 

for  an  isotropic  material,  but  we  assume  that  it  is  equally  valid  for 

propagation  in  an  arbitrary,  anisotropic  material  along  some  specified 

crystal  direction.  The  general  problem  is  then  reduced  to  developing  a 

way  to  calculate  Av  (x.t)  associated  with  that  direction, 
s 

The  surface  wave  penetration  depth  normal  to  the  surface  is 

approximately  the  acoustic  wavelength.  When  the  temperature 

distribution  AT(r,t)  is  essentially  unifomn  over  this  penetration  depth, 

i.  e.  , 9AT/9z  « AT/X  , the  velocity  variation  can  be  approximated  by 

s 

Av  ^ V O'  AT(x,t).  The  velocity  temperature  coefficient  a has 
s so  V V 

been  experimentally  determined  for  a few  common  surface  wave 

materials  including  crystal  quartz  and  lithium  niobate.  Experimental 

I 4 

measurements  of  A(p  for  these  materials  were  compared  with  Eq.  (2) 

using  no  adjustable  parameters,  and  the  quantitative  agreement  was 

within  5%  where  the  theory  was  applicable.  In  those 

experiments,  the  radiative  pulse  widths  were  long  enough  for 

heat  diffusion  to  produce  a uniform  temperature  within  the 

surface  wave  penetration  depth.  In  this  way,  the  need  to  account  for 

propagation  in  the  presence  of  an  axial  temperature  gradient  was 

avoided,  and  the  resulting  analysis  provided  a clear  picture  of  the 

fundamental  surface  wave-thermal  interaction  leading  to  the  induced 
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phase  changes. 

In  general,  depending  on  the  surface  wave  frequency  and  the 
extent  of  thermal  diffusion,  the  temperature  is  nonuniform  over  the 
surface  wave  penetration  depth  as  indicated  in  Fig.  1.  In  the  presence 
of  such  a spatial  inhomogeneity,  surface  wave  propagation  becomes 
dispersive.  The  time  dependent  dispersive  velocity  change  Av^(x,t) 
results  from  a thermal  perturbation  of  the  density  and  the  elastic 
coefficients.  In  order  to  calculate  it  will  be  necessary  first  to  relace 
the  temperature  change  to  the  radiative  absorption  parameters  of  the 
solid  surface,  and  then  to  find  Av^  from  a suitable  dispersion  theory. 

A.  Temperature  Change  Induced  by  Radiative  Absorption 

Radiative  absorption  of  a monochromatic  wave  in  a dielectric 
medium  can  be  described  by  an  absorption  coefficient  0 cm  repre- 
senting the  fractional  intensity  absorbed  per  unit  length,  defined  by 
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B(z')  dz' 


I (r,  z,  t)  = z I {r,  t)  e 
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(3) 


where  I(r,z,t)  is  the  laser  intensity  propagating  in  the  z direction  at 

a depth  z inside  the  surface  of  the  material,  I (r,t)  is  the  intensity  at 

o 

the  surface,  and  z is  a unit  vector  along  the  z-axis.  The  possibility  of 
z dependence  of  the  absorption  coefficient  is  shown  explicitly.  Assuming 
a square  temporal  pulse  of  width  t and  a gaussian  radial  profile  of 
halfwidth  U (HWHM),  we  have 
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0 < t < T 


t > T 


where  P is  the  total  integrated  power  inside  the  material  and  will  be 
related  to  experimentally  determined  quantities  below. 

The  temperature  change  AT{r,t)  resulting  from  the  absorption 
is  found  from  the  heat  equation 
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at 


xV  (AT) 
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(5) 


where  Q = - V I is  the  power  absorbed  per  unit  volume,  p is  the 

density,  is  the  heat  capacity,  and  K is  the  thermal  conductivity  of 

the  material  under  study.  During  the  radiative  pulse  0 < t < t , the 

P 

laser  heating  source  is  given  by 
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For  materials  and  pulsewidths  considered  in  the  present  work,  tempera 
ture  gradients  in  the  z direction  are  much  greater  than  the  radial 
gradients,  with  the  result  that  radial  heat  diffusion  is  insignificant. 
(Solutions  to  the  three  dimensional  heat  equation  based  on  the  work  of 
Olcer  are  consistent  with  this  approximation.)  Consequently,  the 
temperature  change  retains  the  gaussian  radial  dependence  of  the  laser 
profile  and  Eq.  (5)  is  reduced  to  a one  dimensional  problem  by  taking 
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for  which 


AT(r.t)  = ^‘'^T(z.t) 
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where  D = ^’/C  p is  the  thermal  diffusivity  and 
P 


F(z)  = B(z)  exp  I - J B(z')  dz'  | 


A solution  is  obtained  for  an  isolated,  infinite  slab  of  thickness  b with 
the  initial  condition  T(z)  = 0.  Heat  losses  at  the  surfaces  are  neglected, 
yielding  the  boundary  conditions  aT/3z  = 0 at  z = 0,b,  The  formal 
solution  to  Eq.  (7)  is  written  in  the  form  of  a Fourier  series 


T(z,t)  - T (t)  + S T (t)  cos  K z 
o , n n 
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o , n n 

n=l 


where  K = nn/b.  The  equations  for  the  Fourier  coefficients  are 
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The  Fourier  coefficient  F of  the  heating  function  F(z)  depends  on  the 

n 

model  used  to  represent  the  absorptive  properties  of  the  material. 
Several  examples  of  current  interest  will  be  presented  in  Section  III. 


B.  Dispersive  Velocity  Change 

The  velocity  change  will  now  be  related  to  the  temperature 

increase  AT(r,t).  The  starting  point  for  a dispersion  theory  must 

consider  the  spatial  scale  of  the  material  inhomogeneity.  For  exampLO, 

if  the  inhomogeneity  is  concentrated  within  a depth  much  less  than  an 

acoustic  wavelength,  the  boundary  perturbation  can  be  expressed 

8 

concisely  by  an  acoustic  layer  wave  model  . However,  in  our  case 

this  situation  will  only  be  approached  for  the  shortest  times  and  lowest 

irequencies.  For  this  reason,  we  choose  to  represent  the  inhomogeneity 

as  a volume  perturbation  in  the  form  of  a gradual  spatial  variation  of  the 

elastic  properties  induced  by  a temperature  change  for  which 

9 

5^T/3z  < hT/\  . Auld  has  derived  an  expression  relating  the 
^ s 

dispersive  velocity  change  Av  = (v  - v ) to  a volume  perturbation 

s s so 


V CD  O'  ^ 

^ J I -Apu)^  u’^-U  + VgU"  ; AC  : U J dz  (10) 

so  so 


where  u)  is  the  angular  frequency,  U 

vector,  and  V U is  the  strain  field, 
s 

power  flux  in  the  x direction  given  by 


is  the  unperturbed  displacement 

P is  the  unperturbed  acoustic 
s 

OD 


Pj,  j Rl-  Lu  "•  T • x]dz. 
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where  T is  the  stress  tensor.  The  elements  of  the  perturbed 


compliance  tensor  are 


and  the  density  perturbation  is 


(11a) 


hp  = p e AT  ( r , t)  . 
o p 


(11b) 


Here  C°.  . and  p are  the  unperturbed  values,  and  it  is  assumed 
ijkA  o 

e,  e « 1.  In  general,  the  magnitude  of  e will  be  different  for  each 
P 

elastic  coefficient,  but  for  simplicity  the  identifying  subscripts  are  not 
shown. 

In  particular  we  are  interested  in  the  study  of  crystals  exhibiting 
cubic  symmetry  such  as  the  alkali-halides.  The  unperturbed  displace- 
ment vector  for  a wave  travelling  along  the  (100)  direction  of  a non- 

....  10 

piezoelectric  cubic  is  expressed  in  the  form 


where 


U 


u^  (z)  exp  i(u)t  - Kx) 

O 
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and  K = 2ti/X  = (u/ v . The  expressions  defining  the  decay  constants 
s so 
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Of^  and  the  amplitudes  A.  and  are  shown  in  the  Appendix.  Application 
of  the  Cauchy  conditions  readily  extends  the  present  treatment  to 
propagation  on  an  isotropic  elastic  solid.  However,  to  apply  this 
analysis  to  an  arbitrary  propagation  direction  on  an  anisotropic  solid 
of  less  symmetry  would  require  numerical  composite  solutions  for 
the  displacement  vector. 

Using  Eqs.  (7a)  and  (8a)  to  express  AT ( r,  t),  the  following  form 
for  the  velocity  change  is  derived 


Av  (x,t) 
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where  h^,  h^»  and  h^2’  it'  the  Appendix,  are  constants 

independent  of  K involving  the  unperturbed  material  constants,  the 


perturbations  e,  e^,  and  parameters  describing  the  displacement 


amplitudes.  The  factor  u)/4P^  (see  Appendix)  is  also  independent  of  K 


so  that  all  the  dispersive  behavior  is  contained  in  the  lorentzian  terms 


within  the  brackets.  However,  the  dispersive  nature  of  Av  /v  will  be 

s so 


significant  only  if  there  are  terms  for  which  the  Fourier  wave  vector  is 


comparable  to  the  surface  wave  decay  constant,  K =■  o K,  and  for  which 

n i 


the  magnitude  of  is  appreciable.  In  other  words,  the  spatial 
distribution  of  the  teniporaturo  relativt*  to  the  wavt?  penetration  depth 
determines  the  extent  of  dispersive  propagation.  Dispersion  measure- 
ments provide  the  important  experimental  link  to  the  radiative  absorption 
parameters  near  the  surface  because  each  provides  a different  weighted 
sum  of  the  Fitting  the  data  for  several  frequencies  increases  the 

accuracy  with  which  the  and  the  absorption  parameters  are  determined. 

Finally,  we  express  the  induced  phase  change  using  Eqs.  (2)  and 
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in  which  the  gaussian  y dependence  of  AT(r,t)  has  been  used  to  calculate 


f = £ / _IL_ 
B w V in  2 


In  2 ^ 


assuming  a uniform  acoustic  beam  of  width  w.  It  is  interesting 
to  note  that  the  only  dependence  of  the  phase  change  on  the 
laser  halfwidth  CL  is  through  the  error  function  and  this  dependence 


is  exceptionally  weak  for  w/a  >3.  5.  This  form  for  A<|)  is  immediately 


applicable  to  the  general  problem  of  measuring  detailed  surface 
absorptive  properties  by  simply  specifying  F^.  This  will  be  considered 
in  the  next  section. 

III.  DISPERSIVE  PHASE  CHANGE  MEASUREMENTS  OF  8(z) 

In  general,  the  absorptive  properties  of  a solid  within  a depth 

of  several  microns  from  the  surface  may  differ  from  the  relatively 

well  defined  properties  of  the  bulk  material.  This  surface  absorption 

depends  on  surface  preparation,  including  mechanical  polishing  and 

chemical  etching  procedures,  immediate  environmental  impurity 

adsorption,  and  perhaps  on  the  effects  of  surface  dislocations  induced 

during  sample  preparation.  The  resulting  effect  of  these  surface 

perturbations  can  be  described  by  taking  the  radiative  absorption 

coefficient  to  be  a function  of  depth  below  the  surface  3(z).  A simple 

absorption  model  will  be  used  to  represent  the  variation  of  the  absorption 

coefficient;  The  absorbing  surface  region  has  a constant  coefficient 

0 = Bg  cm  ^ which  extends  a characteristic  distance  d into  the  material  at 

which  depth  a transition  to  the  constant  bulK  value  S = 3 cm  ^ occurs, 

B 

This  model  adequately  represents  the  material  absorption  properties 

in  the  region  probed  by  the  surface  wave.  Although  the  sharp  transition 

from  Bg  to  Bg  is  somewhat  unphysical,  thermal  diffusion  during  a laser 

pulsewidth  of  T » 1 usee  causes  the  acoustic  probe  measurement  to  be 
P 

relatively  insensitive  to  the  finer  details  of  the  model. 

Using  this  model  and  Eqs.  (7c)  and  (8b),  it  is  a straightforward 


calculation  to  show 
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Application  to  cases  of  interest  are  effectively  characterized  by 

the  ratio  of  a = S-d  to  8 b which  determines  the  heating  of  the  surface 
S 6 

layer  relative  to  the  bulk  material.  Since  the  penetration  depth  of  the 
surface  wave  is  much  less  than  the  sample  thickness,  the  resulting 
phase  change  is  physically  independent  of  the  parameter  b.  However, 
we  use  the  product  8 b to  define  a dimensionless  quantity  which  can  be 
compared  with  a.  We  emphasize  that  the  variation  of  S^b  should  be 
considered  as  a change  in  the  bulk  coefficient  9^  for  constant  b.  In  this 
section  we  consider  the  dispersive  propagation  resulting  from  two 


examples  of  radiative  absorption;  (a)  a strongly  absorbing  material  for 
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which  3 b » 1,  and  (b)  a highly  transparent  material  having  an 
B 

absorbing  surface  layer.  In  both  cases  experimental  evidence  of  disper- 
sion will  be  presented  and  in  case  (b)  a detailed  study  ot  the  sensitivity 
of  to  absorptive  parameters  is  discussed. 


A.  Strong  Bulk  Absorption;  3.„b  » 1 

Since  the  depth  of  the  surface  layers  under  consideration  are  on 


the  order  of  d ““  1 Um,  the  surface  absorption  a is  weak  even  if  B is 

comparable  to  3 • In  the  limit  of  strong  bulk  absorption  we  have 

B 

3_b  » 1 » a for  which  the  significant  terms  of  the  Fourier  sum  can 
B 

be  approximated  by 
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Thus  I^^  ^Sn  case  and  the  resulting  dispersive  phase  change 

becomes 
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(2o-)^+  (K  /K)^  (a,+  a)^  + (K  /K)^ 

C Vi  i L.  Vi 


(20) 


where  the  radiative  power  inside  the  material  has  been  related  to  the 
Fresnel  reflectivity  R and  the  experimentally  measured  incident  laser 
power  P.  by  P = (l-R)P.. 


The  factors  in  the  sum  given  bv  (3_/K  )^  { 1 - exp  (-K^  Dt)'''/ 

Bn  ri  ' 

2 2 

(B„  + K ) determine  a characteristic  depth  of  the  temperature  distribution, 
5 n 


, which  evolves  in  time  through  thermal  diffusion.  Dispersion  will 

result  when  z is  comparable  to  the  acoustic  penetration  depth  ~ (a  K)~^ 

T i 

Since  z is  initially  « (0„)  • a necessary  condition  for  diserpsive  effects 

T B 

is  ar.Kas  3 . However,  this  is  not  a sufficient  condition  for  all  times, 

1 B 

4 

since  increases  monotonically.  Experimental  evidence  of  dispersion 

fxier  these  Conditions  is  shown  in  Fig.  2 for  surface  wave  propagation 
it  y Tiartz.  Phase  changes  at  acoustic  frequencies  of  12S  MHz 


* r»  in<iu«  I'd  by  the  absorption  of  10.6  am  radiation,  for  which 


0g  241  cm  and  2 and  . 5 for  these  tv^o  frequencies 

respectively.  Although  Eq.  (20)  is  strictly  valid  only  for  a cubic  crystal, 
the  observed  dispersion  on  quartz  shows  qualitatively  similar  effects.  Figure 


shows  the  reduced  phase  change  at  the  time  t per  watt  of  laser  power 

P 

■ ^<l>/(vfgP.)  versus  laser  pulsewidth  The  reduced  phase  change 

is  convenient  in  discussing  dispersion  sinp^  eliminates  the  trivial 

R 

linear  frequency  dependence  in  Eq.  (^5),  and  also  the  uninteresting 

effects  arising  from  different  acoustic  profiles  at  each  frequency.  For 

laser  pulsewidths  ~ 100  msec  the  temperature  distribution  has  become 

essentially  uniform  over  the  31  MHz  wave  penetration  depth  of  ~ 100  um. 

For  these  long  times  the  curves  of  Ab  vs.  t for  124  MHz  and  31  MHz  ■ 

R p 

converge  indicating  insignificant  dispersion.  However,  for  t =2  msec, 

P 

the  temperature  distribution  halfwidth  is  at  a depth  of  only  ~ 80  jjm 

resulting  in  sizeable  dispersion. 

In  this  example,  the  dispersive  phase  change  given  by  Eq.  (20) 

provides  a method  for  measuring  an  exceptionally  large  absorption 

coefficient  without  the  necessity  of  using  an  impractically  thin  sample. 

4 1 

For  example,  an  absorption  coefficient  of  even  3 ~ 10  cm  could  quite 

B 

readily  be  determined  with  surface  waves  in  the  range  of  v 30  - 100  MHz. 


B.  Surface  Layer  Absorption 

Perhaps  a more  demanding  application  is  the  possibility  of 


separating  surface  layer  and  bulk  radiative  absorption  parameters  in 
highly  transparent  materials.  In  particular,  induced  phase  change 


techniques  offer  an  inherently  sensitive  determination  of  the  surface 

absorption  ct  as  will  be  shown  in  the  following  analysis  of  10.  6 ^im 

absorption  on  KCl  surfaces.  The  bulk  absorption  coefficient  at  this 

radiative  wavelength  for  high  purity  KCl  can  range^^  between 

$6  10  to  5 X 10  cm  depending  on  the  growth  procedure.  The 

surface  absorption  can  be  as  low  as  a =“  1 X lO'*^  for  a chemically  etched 

surface  or  as  high  as  a =-  1 - 3 x lO'^  for  a mechanically  polished  surface 

These  ranges  indicate  the  possibility  of  a sample  having  0.1  <■  a/B  b < 10 

B 

and  suggest  that  both  terms  in  given  by  Eq.  (18)  may  be  comparable. 
The  phase  change  is  then  written 


.JfL  ± erf  i p 

V w V za  y V i-R  y t 


^ 1 y i+R 


. / » ^ 1 P,  r "i  , ”2  , ”12  , - -Sb'*’"''! 

' 4P  y b I T 2 2 r 1-e  e 

s (2q'^)  (^“2^  ^“1^  “2^ 


+ E 


n=l  K 


-K^  Dt 

I I„  + I 1 fl  - e ^ 

L Sn  Bn  V y 


I . . "12  . 

^(2trr+(K/K)^  (2«  )^  + (K  /K)^  («,+»,  )^  t (K  /K)^  ' 

1 n 2 n 1 2 n 


In  this  case  the  laser  power  inside  the  material  P = P^(l  + R )/ (I  - R ) has 

/ 

been  expressed  in  .erms  of  P^  the  power  transmitted  through  the  sample 


which  is  a more  appropriate  experimental  quantity  for  a weakly  absorbing 
material.  This  power  term  approximately  includes  the  effect  of  multiple 
reflections.  Further  details  arising  from  the  interference  of  multiply 
[ reflected  optical  field  amplitudes  are  neglected  here  since  the  finesse 

i 

i of  the  sample,  considered  as  an  etalon,  is  extremely  low  for  the 

i ’ 

I 

reflectivity  ~ 4%. 


The  following  analysis  includes  a parameter  study  to  estimate 

the  relative  sensitivity  of  AiJ)  to  a , and  the  individual  parameters  3^ 

and  d as  a function  of  laser  pulsewidth  t . The  logarithmic  derivitives 

d(in  A^)/d(Xn  o)  and  d(Xn  A(J>)/d{^n  3g)  3.re  useful  for  this  purpose  since 

they  directly  indicate  the  fractional  change  in  A<j>  per  fractional  change 

in  a and  B^,  respectively.  In  this  analysis  the  material  properties  of  KCl 

were  used,  for  which  the  acoustic  penetration  depths  (1/e  point)  are 

32  |om  and  128  am  at  the  frequencies  120  MHz  and  30  MHz  respectively. 

The  phase  change  is  evaluated  at  times  t = t for  which  the  dispersion 

P 

theory  is  valid  and  the  derivitives  are  numerically  calculated. 

In  Fig.  3 the  sensitivity  to  surface  absorption  a is  shown  for 
three  values  of  the  ratio  a/B  b = . 1,  1,  10,  keeping  3 = lO’^cm"^  and 

b - , 1 cm  constant.  As  expected,  the  highest  sensitivity  is  approached 
for  shorter  times  and  the  largest  ratio  since  the  heating  is  more 
concentrated  in  the  surface  wave  profile  under  these  conditions.  However 
it  is  important  to  note  that  the  sensitivity  to  a is  quite  high  even  for 
a/Ssb  = 1 at  = 10  msec.  At  this  pulsewidth,  the  temperature  half- 
width is  ~ 200  am  which  is  roughly  a factor  of  6 greater  than  the  surface 


I 


I 


i 


wave  penetration  depth  at  120  MHz.  This  suggests  that  the  experimental 
conditions  under  which  the  integrated  surface  absorption  can  be  readily 
measured  are  well  within  the  technique's  capability.  The  sensitivity  to 
a using  a 120  MHz  surface  wave  is  shown  in  Fig.  3 to  be  only  slightly 
greater  than  for  a 30  MHz  wave  even  though  the  penetration  depth  has 
changed  by  a factor  of  4.  This  indicates  that  the  overall  sensitivity  to 
O is  determined  primarily  by  the  surface  temperature  rather  than  by  the 
temperature  depth  dependence. 

The  sensitivity  to  the  bulk  absorption  coefficient  S is  shown  in 

B 

Fig.  4.  The  logarithmic  derivitive  was  calculated  for  three  values  of 
the  ratio  a/g^b  = . 1,  1,  10,  keeping  a = lo"'^  Og  = 1 cm'\  d = lO"^  cm) 
and  b = . 1 cm  constant.  As  the  pulsewidth  increases,  heat  diffusion 
mixes  the  bulk  and  surface  thermal  energy  in  the  depth  probed  by  the 
acoustic  wave,  resulting  in  an  enhanced  bulk  sensitivity.  It  was 
previously  indicated  in  Sec.  A that  for  strongly  absorbing  materials, 
the  ratio  a/B^b  is  small  and  the  phase  change  essentially  measures 
only  the  bulk  absorption.  The  approach  to  this  limit  is  shown  by 
the  curve  for  a/B^b  = . 1.  Conversely,  when  a dominates,  the  sensitivity 
to  bulk  values  becomes  quite  weak.  It  should  be  emphasized  that  these 
sensitivity  curves  in  Figs.  3 and  4 will  generally  depend  on  the  individual 

magnitudes  of  both  o and  B b and  not  simply  their  ratio. 

B 

Figures  5 and  6 present  phase  change  measurements  obtained  in 
13 

recent  studies  of  surface  absorption  layers  on  chemically  polished  KCl. 
Figure  5 shows  versus  laser  pulsewidth  for  124.  5 MHz  and  28.  0 MHz 


S 1 


Surface  waves  (X^  = 14.  1 i_im  and  62.  5 um,  respectively)  with  theoretical 

fits  of  Eq.  (21)  to  the  data  for  each  frequency.  These  theoretical  fits 

-4 

yield  an  average  value  for  surface  absorption  of  a = (4.  5^  . 5)  x 10  and 

-3  -1 

for  bulk  absorption  0 = (2  ± . 3)  X 10  cm  . The  ratio  a/B_b  = 2 for 

these  data  indicates  a relatively  high  sensitivity  to  a for  short  times  and 

high  acoustic  frequency,  and  low  sensitivity  to  0 except  for  longer  t . 

B p 

Figure  6 shows  the  phase  change  dispersion  in  Ai  versus  laser  pulse- 

£V 

width.  The  data  for  the  two  frequencies  clearly  diverge  from  each  other 

for  short  times,  well  beyond  experimental  uRcertainties  of  about  _+  15%. 

This  divergence  is  direct  evidence  that  absorbing  layers  are  present 

even  after  removing  approximately  10  um  of  surface  material  by  chemical 

polishing.  The  theory  using  the  best  fit  parameters  for  a and  3 repre- 

6 

sents  an  adequate  description  of  the  data  except  for  the  shorter  times 
at  28.  0 MHz.  It  is  believed  this  short  times  discrepancy,  which  is 
emphasized  by  the  log-log  plot,  arises  from  a breakdown  of  the  volume 
perturbation  calculation  of  (Av/v^  ) when  the  characteristic  depth  of 
the  temperature  distribution  becomes  comparable  to  the  acoustic 
penetration  depth  which  is  ~2.2  X^  for  KCl.  Temperature  calculations 
indicate  that  the  ratio  of  z ^ to  the  acoustic  penetration  depth  is  only 
~ 2 at  20  msec  and  decreases  to  ~ . 5 at  2 msec.  Because  of  the  shorter 
penetration  depth  at  124.  5 MHz  (a  factor  of  ~ 4),  the  theory  is  adequate 
at  that  frequency  even  for  the  shortest  times. 

Finally,  it  is  interesting  to  consider  the  conditions  under  which 


the  phase  change  is  sensitive  to  the  individual  parameters  Bg  and  d. 


Thermal  diffusion  is  the  primary  limitation  to  such  a separation  of 


parameters  since  diffusion  eliminates  temperature  differences  over 
distances  ~ d in  a time  d^/D.  In  KCl.  D = . 048  cm^/sec  which 

yields  = 200  nsec,  assuming  a layer  thickness  of  d = lo'"^  cm.  Thus, 
it  would  be  unreasonable  to  expect  the  induced  phase  change  to  be 
independently  sensitive  to  d and  8-  with  laser  pulsewidths  t >1  usee 

^ P ^ 

for  these  conditions.  The  effect  of  thermal  diffusion  is  shown  in  Fig.  7 

which  presents  the  temperature  time  evolution  for  a = lO"  , 

o 

Bg  “ cm  and  b = . 1 cm,  a parameter  set  which  assures  a 

sensitivity  only  to  surface  absorption.  The  temperature  is  calculated 

■ 1 4 1 

for  two  cases:  (i  ) 0g  = 5 cm  , d = 2 X 10  cm  and  (ii)  Bg  = 2 cm"  , 

-4 

d - 5x10  cm,  and  the  percentage  difference  is  shown  as  a function  of 

laser  pulsewidth.  In  Fig.  7a  the  peak  surface  temperatures  are  compared, 

and  in  Fig.  7b  the  depth  of  the  temperature  halfwidth  for  each  case  is 

shown.  Clearly,  for  > . I msec  the  temperature  distribution  depends 
\ 

only  on  which  is  the  same  for  both  cases.  However,  a Q-switched 
€©2  laser  with  . 5 ijsec,  or  a TEA  laser  with  t . 1 tisec  (total 

energy  width),  would  certainly  be  useful  to  observe  the  individual 
dependence  on  the  parameters  0 and  d in  KCl  surface  layers. 


IV.  CONCLUSION 

The  present  theory  describes  the  dispersive  behavior  of  the 
thermal-acoustic  interaction  and  properly  accounts  for  this  effect  on 
the  magnitude  of  the  laser-induced  phase  change  of  a surface  acoustic 

.S5 


wave.  Although  the  analysis  here  has  stressed  cubic  materials,  it  is 
quite  straightforward  to  extend  the  present  treatment  to  materials  of 


arbitrary  symmetry.  In  addition,  the  possibility  of  including  a material 
inhomogeneity  on  a scale  much  less  than  the  acoustic  wavelength  requires 
only  minor  revision.  Experimental  measurement  of  dispersive  surface 
wave  propagation  offers  the  possibility  of  examining  the  detailed  charac- 
teristics of  radiative  absorption  at  material  surfaces.  The  facility  of 
these  techniques  is  enhanced  by  the  time  evolving  thermal  gradient 
which  adds  an  experimental  degree  of  freedom  reducing  the  need  for 
measurements  at  many  different  wavelengths.  In  addition,  the  experi- 
mental sensitivity  inherent  in  phase  measurements  increases  the  utility 
of  this  technique. 

We  have  considered  two  examples  for  analysis  here:  strong  bulk 

absorption,  and  surface  layer  absorption  in  essentially  transparent 
materials.  The  latter,  a current  problem  in  high  power  laser  window 
technology,  was  examined  more  thoroughly  to  emphasize  the  sensitivity 
of  phase  change  measurements  to  the  intricate  details  of  absorbing 
surface  layers.  Experimental  results  were  presented  indicating 
dispersive  propagation  of  surface  waves  on  KCl  following  absorption  of 
10.6  um  radiation.  Fitting  the  theory  to  the  data  provided  a measure- 
ment of  the  integrated  surface  absorption  o as  well  as  the  bulk  absorpLion 
coefficient  0 cm  ^ . 

D 

The  application  of  this  measurement  technique  to  examine  a wider 
range  of  surface  piienv>niena  may  be  i unside  red.  Certiiinly  the  layer 


r ' ^ 

absorption  model  is  particularly  appropriate  to  describe  absorption  in 
a thin  dielectric  or  metallic  film  deposited  on  an  otherwise  transparent 
substrate.  This  relates  to  the  need  to  measure  a for  highly  reflecting 
films  and  protective  film  coatings  associated  with  high  power  laser 
optics.  The  possibility  of  obtaining  the  spatial  distribution  of  an 
absorption  coefficient  3(z)  near  a surface  may  be  particularly  interesting 
if  the  absorption  is  directly  related  to  the  diffusion  of  impurities  into  the 
surface.  In  this  case  a tuneable  laser  resonant  at  the  impurity  absorption 
would  induced  a thermal  gradient  sensitive  to  the  particle  diffusion 

properties.  The  availability  of  tuneable  infrared  lasers  also  suggests  j 

the  possibility  of  tuning  through  restrahl  resonances  at  the  surface  of  a S 

! 

solid.  Surface  waves  propagating  on  the  solid  surface  would  serve  to  | 

I 

! 

detect  these  very  high  crystal  absorptions  without  the  need  for  fabricating 
thin  samples.  The  theory  developed  in  this  paper  establishes  a starting 
point  for  the  extension  of  these  surface  wave  techniques  to  the  study  of  a 
variety  of  surface  processes  with  increased  precision  and  sensitivity. 

The  authors  wish  to  thank  S,  Allen  and  A.  Timber  of  Hughes 

i 

Research  Laboratories  for  the  KCl  samples  and  their  preparation.  They  j 

1 
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fabrication,  and  D.  Smith  for  his  invaluable  assistance  in  writing  the  | 

computer  programs. 
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APPENDIX:  DISPERSIVE  PROPAGATION  IN  CUBIC  MATERIALS 


It  is  the  purpose  of  this  Appendix  to  derive  an  expression  for  the 
change  in  surface  wave  velocity  induced  by  a subsurface  thermal  gradient 
in  materials  having  cubic  symmetry.  The  unperturbed  solution^^  for  the 
mechanical  displacement  vector  will  be  presented.  The  velocity  change 
can  then  be  expressed  in  terms  of  this  solution  and  the  existing  tempera- 
ture distribution. 

We  begin  with  the  wave  equation  for  the  displacement  component 
of  an  elastic,  anisotropic  medium,  in  the  absence  of  external  body 
forces  and  piezoelectric  effects: 


''ijki 


9x . dx . 3 X , 

J I I 


(i,j,k,f  = 1,2,3)  (Al) 


The  presence  of  the  second  term  accounts  for  possible  spatial  variation 
in  the  elastic  constants.  In  the  present  work,  this  variation  is  expressed 
in  the  form 


'ijk.4 


'ijki 


1 + 


ijk;. 


t) 


(A2) 


and  the  density  variation 

p = p jl+e  AT(r,t)"i  . 
o L P 


( 

Here  and  are  tlie  unperturbed  values,  and  the  dependence  on  the 

—A 

temperature  increase  AT(r,t)  is  weak  since  the  coefficients  e , 

P 

,4 

e...  . 10  /°K  for  materials  under  consideration.  We  consider 

ijkjt  ~ 

surface  wave  propagation  along  a crystal  symmetry  direction  taken  as 

the  X axis,  and  the  normal  to  the  surface  is  parallel  to  the  z axis.  We 

assume  the  x and  z variation  of  AT(r,  t)  is  slow  enough  that 

3c.  /3x  « 2nc°,  /\  and  3c.  /3z  « 2nc°,  /X  , where  X is 

ijki  ijk.i  s ijk£  ijkjt  s s 

the  acoustic  wavelength.  In  this  case  the  second  term  in  Eq.  (Al)  may  be 
neglected.  This  condition  is  readily  satisfied  in  the  presence  of  such 
small  temperature  coefficients.  We  have  then,  for  slowly  varying 
elastic  properties 


P 


5^, 


'ijkjt 


3x.  3x 
J I 


(A4) 


Solutions  which  represent  a wave  propagating  parallel  to  the  x-axis, 
decaying  with  increasing  distance  z from  the  surface  z = 0 and  having  no 
dependence  on  the  y-coordinate  are  expressed  in  the  form 


U 


u^  (z)  exp  i(u)t  - Kx) 

O 

-i  u^(z)  exp  i (cut  - Kx) 


(A  5) 


whe  re 


- Kz  - Qf^  Kz 

+ A^e  ' 


■(A  6) 


Ul(z)  = A^e 


-a 

u^(z)  - B^c  r B^e 


-Of,  Kz 


for  a surface  wave  of  angular  frequency  u) , unperturbed  velocity 

and  wave  vector  K - 2n/X  = uj/ v . The  decay  coefficients  are  determined 

s so 

by  the  pair  of  equations 


af+  = - j (1  - Y R)  + Y^(l  - R)  - (Y  + y')^  ] 

1 2 Y L J 


(ot^a^)  = (1  - Y R)  (1  - H) 


where  Y = /c°  , y'  = /c°  and  R = p /c°  and  the  unperturbed 

44  1 1 12  11  o so  44 

o 

elastic  constants  are  expressed  in  reduced  subscript  notation  The  three 


amplitude  ratios 


h = a:  = 

I 


a . (y  + Y ' ) 
Y (1  - R)  - 


(i  = 1,2) 


(A  10) 


“2  ^^2 


(All) 


specify  the  solution  to  within  a constant  determined  by  the  total  mechnical 
energy.  Finally,  the  unperturbed  acoustic  power  flux,  , is  given  by 


® - ^ro  o o n 

- — I c , + c b (Q'+b  )-c  Qr  b I 

u)  a^Lll  44  11  1 12  11_ 
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o,  L 11  44  2'  2 2'  12  2 2 _ 


I r,  _o  ^ 


“l*“2 


^'ll  '■'=44['=l'“2+‘’2'  +''2'°'i+N>] 
-'=12'“l‘’l*“2‘’2']- 


(A  12) 


The  dispersive  velocity  change  induced  by  the  presence  of 
AT(r,t)  is  given  by  Eqs.  (10)  and  (14)  under  the  conditions  discussed  in 
Section  IIB.  For  the  surface  wave  solution  given  above,  the  constants 
hj^,  h^  and  h^^  given  by 


i = 2a  j d - v^  p c (1  + bf ) 1 
1 iLl  so  op  I J 


(A  13) 


h-  = 2a  a^  j d - V ^ p e (1  +b^)"', 
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h = 2(a  + o ) a j d - v ^ p e (1  + b b ) 1 
12  121L12soop  12J 
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where 
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11*11  J ■ ‘'44"44<^2'  ^2^' 


'^12~  ^11*11  [^^“l^l^z'^z]  ■''12*12  “2^2 


(A  18) 


This  completes  the  expressions  needed  to  calculate  the  dispersive 
velocity  change  ior  wave  propagation  along  a symmetry  axis  in  cubic 
materials.  Equations  (A  13)  - (A  18)  and  Eq.  (21)  have  been  applied  to 
calculate  phase  changes  induced  by  radiative  absorption  at  KCl  surfaces. 
The  results  of  these  calculations  are  shown  in  Figs.  5 and  6. 
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FIGURE  CAPTIONS 


Fig.  1 General  experimental  configuration  for  the  study  of 

radiatively  induced,  acoustic  surface  wave  velocity 

dispersion.  Incident  laser  radiation  is  absorbed  in  the 

shaded  region  which  results  in  a temperature  change  having 

a depth  distribution  indicated  by  AT(z,t).  This  temperature 

distribution  is  probed  by  an  acoustic  surface  wave  of  wave 

vector  K , whose  energy  distribution  is  shown  having  a 
s 

penetration  depth  comparable  to  the  acoustic  wave  lengtn 

X . Coordinate  axes  shown  are  not  intended  to  coincide 
s 

with  the  crystal  axes  of  the  sample. 

Fig.  2 Phase  change  dispersion  on  Y-cut  quartz  as  a function  of 

laser  pulsewidth.  The  reduced  phase  change  per  watt  of 

incident  laser  power,  = A(^/(vf„P. ),  is  shown  for 

R B 1 

acoustic  frequencies  of  3 1 MHz  and  124  MHz.  Dashed 
curves  are  shown  as  guides  to  the  eye. 

Fig.  3 The  relative  sensitivity  of  the  induced  phase  change  to  the 

integrated  surface  absorption  o is  indicated  by  the  logarithmic 

derivitive  d(in  A(^)/d(in  a)  as  a function  of  laser  pulsewidth. 

This  is  shown  for  three  values  of  the  parameter  ratio 

-3  -1 

(j/Sgb,  keeping  = '0  cm  and  b = . 1 cm  constant. 

Curves  are  shown  for  acoustic  frequencies  of  120  MHz 
and  30  MHz. 

Fig.  4 The  relative  sensitivity  of  the  induced  phase  change  to  the 

bulk  absorption  coefficient  8 is  indicated  by  the  logarithmic 
derivative  d(£n  A^/d(/n8g)  as  a function  of  laser  pulsewidth. 
This  is  shown  for  three  values  of  the  parameter  ratio 
a/8gb,  keeping  a = lO"'^  (8g  = 1 cm'\  d = lo"^  cm)  and 
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b = . 1 cm  cmstaxit.  Curves  are  shown  for  acoustic 
frequencies  of  120  MHz  and  30  MHz. 

Fig.  5 Phase  change  data  versus  laser  pulsewidth  for  a 

chemically  polished  KCl  surface  at  acoustic  wave 
frequencies  of  124.  5 MHz  and  28.  0 MHz.  Bars  indicate 
± 15%  experimental  uncertainty.  Solid  curves  are 
theoretical  fits  to  the  data  for  each  frequency.  The 
average  value  of  a and  8 obtained  from  the  best  fit  are 
indicated. 


Fig.  6 Phase  change  dispersion  versus  laser  pulsewidth  is 

indicated  by  the  increasing  divergence  of 

decreases.  These  data  were  taken  on  a chemically  polished 
KCl  surface  at  acoustic  wave  frequencies  of  124.  5 MHz 


and  28.  0 MHz;  experimental  uncertainty  of  ± 15%  is 

indicated  for  each  data  point  by  the  bar.  Solid  curves 

are  theoretical  fits  to  the  data  for  both  frequencies.  The 

average  value  of  a and  obtained  from  the  best  fit  are 

B 

indicated. 

Fig.  7 Measurement  sensitivity  to  individual  values  of  absorption 

parameters  8g  and  d depends  strongly  on  thermal  diffusion. 
This  is  indicated  in  (a)  by  comparing  the  peak  surface 

temperatures  versus  laser  pulsewidth.  T_  is  calculated 

-1  -4  ^ -1 

using  83  = 2 cm  , d = 5 x 10  cm  and  T^  using  83=  5 cm  , 

d = 2 X 10'^  cm.  In  (b)  the  halfwiddi  of  the  temperature 

depth  distribution  versus  laser  pulsewidth  is  shown  for 

these  same  parameter  sets. 
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